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f^ Abstract: We investigate how relativistic, nonabelian plasmas approach equilibrium in a 

i ■ general context. Our treatment is entirely parametric and for small Yang-Mills coupling a. 

First we study isotropic systems with an initially nonequilibrium momentum distribution. 

^G . We consider both the case of initially very high occupancy and initially very low occupancy. 

Then we consider systems which are anisotropic. We consider both weak anisotropy and 
large anisotropy, and allow the occupancy to be parametrically large or small. Writing 
the typical momentum of an initial excitation as Q and the final temperature as T nna i, full 

(— «) equilibration occurs in a time t eq ~ a_2 ^fin a l f° r tfinai > Q, and t cq ~ a_2( 3 5 ^final ^ or 

^O | 7fi na i < Q, unless the initial system is sufficiently anisotropic and T nna i > a^Q, in which 

o 
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case equilibration occurs somewhat faster, t eq ~ max(a T , a t Q i T^J,). 
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1 Introduction and summary of results 

Nonabelian plasmas out of equilibrium are a rather generic feature of early universe cos- 
mology. For instance, they can arise in the process of reheating after inflation, whether by 
perturbative decay of inflatons [1], by resonant decay of an inflaton into Standard Model 



fields [2], or due to other late-decaying relics. In some cases the nonabelian plasmas created 
in these processes can be locally very anisotropic; for instance, Asaka et al [3] show that 
late decays of inflatons can produce jet-shaped regions of heated plasma which may be 
important for electroweak baryogenesis; the interiors of these jets are generically highly 
anisotropic. Nonequilibrium and anisotropic plasmas can also arise in cosmological phase 
transitions, such as the electroweak phase transition [4]. In all of these cases, due to 
asymptotic freedom in QCD and the weakly coupled nature of weak isospin, the relevant 
nonabelian coupling is relatively small. 

The early stages of highly relativistic heavy ion collisions also probably produce a 
nonabelian (QCD) plasma out of equilibrium. Probably in real-world applications the 
coupling is not small and a perturbative approach is suspect; but at least in the theoretically 
clean limit of extremely high energy collisions the coupling is also weak in this context. 

In any case we think it is a well motivated theoretical question to ask, how in general 
does a nonabelian plasma approach equilibrium? For the case of a plasma which is relatively 
close to equilibrium, we believe that this is well understood. The dominant processes in this 
case are elastic scattering and inelastic, number-changing (splitting and joining) processes. 
These are well described by an effective kinetic theory [5]. There is a natural energy scale 
T, and the time for the system to approach equilibrium is parametrically t cq ~ a~ 2 T~ l up 
to logarithmic factors which we will systematically ignore. 

However, it is less clear what physics is most important for a system which starts out 
very far from equilibrium. This is particularly the case if the system is also anisotropic. 
In this case, there can be plasma instabilities [6-11] which may dominate the dynamics in 
some cases. To our knowledge, there has been no comprehensive study of this case. The 
literature is also somewhat fragmentary for the case of a system which is isotropic, but 
in which the modes which dominate the system's energy have typical occupancies much 
larger or much smaller than 1. 

In this paper we will give a parametric treatment of each of these cases; systems 
which are far from equilibrium but isotropic, and systems which are far from isotropic. 
By a parametric treatment we mean that we will identify the most important physics 
in each case, as determined by counting powers of the gauge coupling a, and estimate 
relevant momentum and time scales in the process of equilibration as powers of a. We 
will not attempt to keep track of factors of the logarithm of the coupling. Nor will we 
attempt to estimate any order-1 coefficients, or determine the range of a values over which 
our parametric estimates are reliable. If a more quantitative estimate is needed in any 
particular situation, then we have at least identified what the relevant physics is in making 
such an estimate. 

We will not present estimates of the equilibration of systems which are initially spa- 
tially nonuniform (inhomogeneous) . This presents too general a class of situations to allow 
any comprehensive study. However there is a general strategy for using our results to 
study these systems as well. Namely, one can consider the way in which spatial nonuni- 
formity leads to local momentum-space nonuniformity through propagation, and then use 
our results for the evolution of momentum-space anisotropic systems to determine the local 
physics which this causes. We will use the results of this paper to study one particularly 



interesting case, that of a system under 1-dimensional "Bjorken" expansion, in a future 
publication. 1 

In the main body of the text we will go through each case in detail, introducing the 
physics relevant to each situation as the need arises. However, for the convenience of the 
reader we will begin with a summary of our results. 

First we consider isotropic plasmas. For simplicity we only consider plasmas where all 
physics is initially dominated by excitations which have a single characteristic momentum 
scale Q, with a typical occupancy f(Q) which we write parametrically as f(Q) ~ a~ c . 
Cases with c > represent initially overoccupied systems; those with c < represent 
initially underoccupied systems. 

We find three cases. For c > 1 (extreme overoccupancy) the Nielsen-Olesen instability 
[13] rapidly converts the system into one with a larger Q and c = 1. 

For < c < 1 (overoccupancy), the dominant processes are elastic scattering and 
number-changing effective 2 — > 1 "merging" processes. Excitations quickly arrange into an 
f(p) oc T*/p distribution with a cutoff scale p max . Initially p max ~ Q arid T* ~ a~ c Q, but 
both scales evolve with time, see eq. (2.30): 

2-2c „ 8 1 -6-c „ 4 3 , oio„„ix , 

Pm^-a — Qrtr , T*~a—Qrt-7 , (tXx-^Q- 1 ) . (1.1) 

Equilibration occurs when p max ~ T* ~ a~ c > 4 Q, at £ C q ~ a~ 2+ iQ~ 1 ~ a~ 2 T _1 . The 
final temperature is determined by the initial energy density; and t eq is the characteristic 
equilibration time for a bath of temperature T. 

For c < (underoccupancy) , the physics is a little more involved (though a basically 
correct exposition was given in Ref. [12], in a slightly different context). The most im- 
portant physics is the formation of a bath of low-momentum "daughter" particles, which 
thermalize and begin to dominate the system's dynamics at a time scale t ~ a '"sQ -1 . 
The soft thermal bath then catalyzes the breakup of the hard excitations, absorbing their 
energy into the thermal bath. Both processes are dominated by number-changing effective 
1—7-2 processes, and it is essential to include the modification of the number-changing 
rate due to the Landau- Pomeranchuk-Migdal (LPM) effect [14, 15]. The temperature of 
the thermal bath rises with time as T ~ a i ~ c Q^t 2 . The hard excitations are consumed and 



thermalization completes when T = Tg na i ~ a * Q at t cq ~ a 2+ s Q 1 . The temperature 



3c 
■-eq 

is again dictated by the initial energy density; the equilibration time is t C q ~ a a_2 ^finai' 
longer by a c ' 8 than the characteristic equilibration time of a thermal bath at temperature 
Tfinal- But t eq is much shorter than the large- angle elastic scattering time scale for the 
initial distribution of excitations. Therefore, large angle change through elastic scattering 
between the initially present excitations actually plays no role in the system's thermaliza- 
tion. 

There is a simple way to understand the equilibration time for an underoccupied 
system. Because of LPM modified bremsstrahlung emission, a hard p S> T excitation 



1 We cannot help giving away one central result of our future publication. Assuming the energy density 
scales as e ~ aT 1 Q i r~ 1 before equilibration, and using the estimate for the final equilibration time given 
in the abstract, one finds r cq ~ a^^Q -1 , slightly different than the conclusions of Baier et. al. [12]. 



Region 


Boundaries 


timescale for 
new physics 


New physics 
which occurs 


1 


c> 0, 3d - c + 1 > 0, 7d - c + 1< 


-d-Y + Yc . 


joining reduces anisotropy 


2 


c<0, 3d + c+l>0, 5d + c+l<0 


-7d-7-3c -, 

a 4 Q 


new plasma instabilities 


3a 


c < 0, 5d + c + 1 > 0, d < 0, Id + 3c + 1< 


-lld-5-3c „ -, 

a 2 Q 1 


new plasma instabilities 


3b 


d > 0, 8d + 3c + 3 > 0, 9d + 3c + 1< 


-35d-lB-9c -, 

a 6 Q 


new plasma instabilities 


4a 


7d - c + 1 > 0, Id + 3c + 1 > 0, d < 


3d-3+3c _ -, 

a 2 Q 1 


angle randomization 


4b 


d > 0, 23d + 9c + 3 > 0, 3d + c - 1< 


5d— a+ac -| 
a 2 Q 


angle broadening 


5 


3d + c + 1 > 0, 8d + 3c + 3 < 


a 2 Q x 


new plasma instabilities? 


6 


9d + 3c + 1 > 0, 23d + 9c + 3 < 0, 9d + 3c - 1< 


-4d-6 i 

a 3 Q i 


new, weak instabilities 


7 


3d + c-l>0, d-c+l>0 


a 2 Q x 


angle broadening 


8 


c> 1, d — c+1 < 


less than Q _1 


Nielsen-Olesen instabilities 


9 


3d + c + 1< 0, c < or 3d - c + 1< 0, 0<c<l 


as per isotropic 


as per isotropic 


10 


9d + 3c - 1 > 0, 37d + 15c + 7 < 


-7d-'25+lic , 

a 12 Q i 


thermal bath elastic scatt 



Table 1. Main regions in the occupancy (c)— anisotropy(d) plane, indicating the time scale on which 
the physics changes and the nature of the new physics. 



in a thermal bath loses energy at a rate dE/dt ~ a 2 T 2 y / E/T [15]. Therefore £ e q is 
the characteristic time scale for an excitation of energy Q, traversing a thermal bath of 
temperature Tfi na \, to radiate away its energy. Equilibration could not proceed faster than 
this. 

Having treated isotropic systems, we then turn to anisotropic ones. Again we assume 
that the system starts with excitations possessing a single characteristic momentum scale 
Q. But the distribution of particles is now characterized by two quantities; a typical 
occupancy, f(p) ~ a~ c , and a measure of anisotropy characterized by a strength a~ . For 
weakly anisotropic systems with d < 0, we define e = a~ as the relative variation of the 
occupancy with direction (so if e = 0.1 then the occupancy in some directions is 10% larger 
than in other directions). For strongly anisotropic systems with d > 0, we define 5 = a~ d , 
as the angular range within which most of the excitations reside. We concentrate on the 
case of an oblate distribution, with most excitations' momenta p lying within an angle 5 
of the xy plane. In this case, if 5 = 0.1 then the occupancy is f(p) ~ oT c if \p z \ < 0.1|p|, 
but f(p) is small outside this range. We also consider prolate distributions with most 
excitations' momenta within an angle 5 of the z axis. 

Our findings for oblate anisotropic plasmas are summarized in figure 1 and table 1. 
Depending on the values of c and d (occupancy and anisotropy), the system's evolution 
can take many different forms. In almost all regions plasma instabilities play a key role in 
the dynamics; but after some time scale new physics comes to dominate and the behavior 
of the system changes. We have only tried to track the evolution until the first significant 
change in the dominant physics. In summary, the regions are: 

1. Plasma instabilities induce joining processes which suppress the anisotropy. 



Anisotropy 

A 



d=ln(8)/ln(<x) 



: Broadening: non-HTL, 
unsaturated plasma instabilities 




9: Plasma instabilities 
Play No Role 



Occupancy 
c=-ln(f)/ln(a) 



8:Nielsen- 
Olesen 
Instabilities 



Figure 1. Main regions in the occupancy-anisotropy plane. Descriptions of each region are in the 
text, and in Fig. 13. 



2. Plasma instabilities induce splitting processes, split daughters are anisotropic and 
generate new plasma instabilities which come to dominate the dynamics. The most 
important daughters have occupancies limited by re-absorption. 

3. Same as 2. except re-absorption is irrelevant; the most important momentum scale is 
set by the angle randomization of split daughters. 

4. Plasma instabilities randomize the directions of hard excitations. 

5. Radiated daughters induce new plasma instabilities before the original plasma insta- 
bilities finish growing to their "saturation" size. This region is not fully understood. 

6. Radiated daughters form a nearly-thermal bath, but residual anisotropy of the bath 
generates new plasma instabilities. 

7. Plasma instabilities broaden the momentum distribution before the unstable modes' 
growth can saturate. The hard-loop approximation does not apply in this region. 

8. Nielsen-Olesen instabilities rapidly lower the maximum occupancy. 

9. Plasma instabilities play no important role. Equilibration proceeds as in an isotropic 
system. 

10. Plasma instabilities cause splitting which generates a soft thermal bath. The bath 
comes to dominate the dynamics through elastic scattering. (This region lies off the 
edge of Fig. 1.) 

Note that in regions 1, 4, 7, and 8 the processes which redistribute particle momentum (q) 
become less efficient with time; whereas on the contrary, in regions 2, 3, 5, 6, and 10 the 
new physics which emerges actually makes the redistribution of particle momentum more 
efficient with time. 



Though we have not followed the evolution of anisotropic systems through to final 
thermalization, we have an estimate for the final thermalization time. First one determines 
Tfinal ~ e 3 ; for d > it is T fina i ~ a^rQ, for d < it is T final ~ a~fQ. If T final > Q then 
t cq ~ a- 2 T^ v If Tflnai <C Q then t cq ~ cTTQ? T fl "^ if T final > a§Q and t cq ~ oT 2 Q^T^ 

2 

if Tfinai < ct^Q. Both estimates are the time it takes for an excitation of momentum Q 
to lose its energy while traversing a thermal bath of temperature ~ T nna i; the difference is 

2 

that for T nna i > a 3 Q, the thermal bath is somewhat anisotropic and plasma instabilities, 
not elastic scattering, determine how the p ~ Q momentum excitation loses energy. If the 
resulting i eq is shorter than ci~ 2 T^ nsl , then the hard p ~ Q excitations break up before the 
remaining thermal bath can equilibrate with itself. But it still takes time ~ a_2 tfinai ^ or ^ e 
bath to fully self-equilibrate, so this is remains a lower bound on the total thermalization 
time. 

Having summarized our findings, we will now present the details leading to these 
results. Since we presented the main conclusions here, we will not end with a discussion 
or conclusions section. 

2 Isotropic distributions 

In this section we will consider systems which are both homogeneous and isotropic. Never- 
theless they can be very far from equilibrium. Such systems might arise rather generically 
during preheating after inflation. We will also re-encounter many of the physical pro- 
cesses relevant to this case when we consider systems which are not isotropic. We begin 
by reviewing the relevant physical processes, then we will consider their application to 
nonequilibrium systems. 

2.1 Elastic and inelastic scattering 

We begin by reviewing the basic physical processes relevant for equilibration in a plasma. 
All of the physics discussed here is well known, and is presented in more complete detail 
in Ref. [5]. At weak coupling and at sufficiently long length and time scales, the dynamics 
of the plasma can be described by that of long-lived weakly interacting quasiparticle exci- 
tations. The perpetual interaction with the medium affects the quasiparticle's dispersion 
relation, and to leading order, the correction looks like an effective mass 

E(p) ~ yV + ™ 2 , (2.1) 

which introduces a new scale, the screening scale m, which is parametrically 

2 _ r hp) 



m z ~ q / z-2-1-, (2.2) 

J P P 

where J" is a shorthand for J d 3 p. 

Under certain conditions, the non-equilibrium dynamics of the plasma can be described 
by the Boltzmann equations that schematically read 

(d t + v(p)-V x )f(x,p,t) = -C[f], (2.3) 






Figure 2. Diagrams relevant in elastic scattering of gluons. 

where f(x,p,t) is the phase space density of quasiparticles in the plasma, and v(p) is the 
velocity of a quasiparticle with momentum p. In the following, we will restrict ourselves 
to spatially homogeneous systems so that / is only a function of p and t. In the presence 
of several particle species (gluons, up- and down-quarks, etc.) the distribution function is 
a multicomponent vector with a component for each species. In the absence of external 
forces, the quasiparticles change their momentum states by mutual interactions; C[f] is a 
spatially local collision term that represents the rate at which particles get scattered out 
of state p minus the rate they get scattered into this state. 

The conditions under which the non-equilibrium dynamics of a theory can be described 
by the Boltzmann equations of eq. (2.3) are 

• that the typical size of the wave packets of quasiparticles is smaller than the mean 
free path of quasiparticles, and 

• the quantum mechanical formation time of scatterings is small compared to the mean 
free time. 

Generally these conditions will be met for those degrees of freedom in a plasma with typical 
occupancy f(p) <C a -1 . An exception is that certain inelastic processes can have formation 
times which exceed the mean free time between scatterings for the excitations involved. 
This case requires special treatment (the LPM effect), which we will return to at length 
below. 

There are two qualitatively different processes that contribute in leading order to the 
collision term, elastic scattering and near-collinear splitting processes. The collision term 
for elastic scattering reads 

C 2 ^[f]= \j \M(p,k,p>,k')\ 2 (2ir) 4 5^(P + K-P'-K') (2.4) 

z Jk,p',k' 

x{f(p)f(k)[l±f(p')][l±f(k')] - f(p')f(k')[l±f(p)][l±f(k)]}. (2.5) 

Here P, P' , K, and K' denote on-shell four-vectors. The [l±/]-factors arise from final state 
Bose stimulation or Pauli blocking. M is the elastic scattering amplitude in non-relativistic 
normalization, related to the usual relativistically normalized matrix element A4 by 

{2p ){2k ){2p' ){2k' ) 



In the leading order for gauge bosons, the M. originates from diagrams in Fig. 2 (from 
now on we will only consider gauge bosons, which due to the stimulated interaction rates 
dominate the equilibration). 2 In vacuum the matrix element reads 

, yl9 o / su st tu\ , ... 

N (vacuum ~ <** [3 ~ ~p ~ ^ ~ j J , (2-7) 

where s, t, and u are the usual Mandelstam variables. The process is dominated by small 
momentum exchange, and in this limit the matrix element becomes 

2 

IMI 2 (2 8) 

l 1 '- 1 Ivacuum /■ 2 \2 ' V^ ,u 7 

where q± is the momentum transfer in the elastic scattering. In vacuum the total cross- 
section ~ f d 2 q±\M\ 2 is infrared divergent, but this divergence is regulated by including 
the medium dependent self-energy to the exchange line in the diagram of Fig. 2. Then the 
amplitude becomes [16] 3 

2 
l M ' 2 ~ 2,2°, 2V ( 2 ' 9 ) 

A particle traveling through the plasma undergoes successive uncorrelated elastics 
scatterings that diffuse the momentum of the traveling particle. The momentum of the 
particle evolves as a random walk in momentum space, characterized by the momentum 
diffusion constant (/elastic! the average squared momentum transfer Ap 2 grows linearly in 
time, 

A/ ~ ^clastic*- (2.10) 



The rate of elastic scatterings is 

aT P i 



d 2 qi 



dq z f \M\ 2 f(p')[l+f(p'-q)], (2.11) 

Jp' 



which for qj_ < \p\ is 



dT e \ a 



2 



d 2 q± q\{q\ + 



m' 



f(p')[l+f(p')}, (2-12) 



p 



and the mean squared momentum transfer per unit time is 



r/r c i 



Elastic rJ I d g ± -2 — q ± . (2.13) 



2 If the typical initial occupancy is large, we must be speaking of gauge bosons. If the initial occupancy 
is small and the system is composed primarily of fermions, then the radiated "daughter" particles are 
primarily gauge bosons. The rate of radiation differs only by an order- 1 group Casimir factor and all 
parametric estimates are identical to the case with primary gauge bosons. All of these remarks also apply 
to the anisotropic case. 

3 The screening does not affect nearly static magnetic fields, which is why the matrix element scales as 
l/<7x for very small q±. Because of this behavior the total cross-section actually remains logarithmically 
divergent, but the divergence is too weak to influence the rate of angle or particle number change. 




Figure 3. Diagrams responsible for inelastic scattering processes in the plasma. The thick lines 
represent an initial particle and its two nearly collinear daughters; the thin lines show a scatter- 
ing process which introduces a small amount of momentum, rendering the splitting kinematically 
allowed. 



Collisions with small q± ~ m are the most frequent but change momentum of the propa- 
gating particle the least, while scatterings with large q± are rare but produce larger angle 
deflections. Inserting eq. (2.12) into eq. (2.13), one sees that all available scales of momen- 
tum transfer contribute equally to the momentum diffusion, such that ^elastic reads 4 



Elastic ~ a 2 / /(p)[l+/(p)]. 



(2.14) 



Due to the peculiarities of soft and collinear enhancements in QCD, there is another 
class of processes which are also leading-order. These are the inelastic, number chang- 
ing, collisions [5, 17]. In the vacuum, 1 -H- 2 processes are kinematically disallowed for 
massless particles, but in the presence of the medium, soft elastic scatterings can take the 
particles slightly ~ m off-shell allowing for subsequent "1 -f-> 2" nearly collinear splitting 
processes, depicted in Fig. 3. A particle propagating through a plasma experiences soft 
elastic scatterings with rate 



rei~^/ f{ P ')[l+f{ P ')] ~ ^ 



(2.15) 



In vacuum, each individual uncorrelated scattering event induces a collinear radiation with 
probability a per logarithmic range in angle, per logarithmic range in energy. The same 
applies in-medium (modulo Bose stimulation factors) provided that the formation time 
of the emission process is shorter than the mean time scale between scattering events. 
When the formation time is longer than the time between scatterings, there is destructive 
interference between emission processes, the LPM effect, which we will discuss in due 
course. Assuming the formation time is shorter, and neglecting the log over emission 
angles, the particle splitting rate is parametrically 



UL split 

dp/p 



«r e i, 



(2.16) 



where p is the energy of the emitted particle. (The superscript BH stands for Bethe- 
Heitler, since the regime where coherence effects are negligible corresponds to the original 



That all the available momentum transfer scales contribute to the momentum diffusion gives rise to a 
logarithmic enhancement, which is however neglected in our power counting. 



calculation of Bremsstrahlung radiation in QED by Bethe and Heitler [18].) Including the 
statistical factors, the effective collision term for "1 <R- 2" processes reads 

r dT BH 
C" 1 " 2 " [/] ~ J dk^± (/(p) [l+f(p-k)] [l+f(k)} - f(p-k)f(k) [l+/(p)]) , (2.17) 

where energy and momentum conservation are implied. 

2.2 < c < 1: Cascade by elastic scattering 

We start by considering isotropic and homogeneous systems whose initial distribution at 
t = has occupancies ~ aT c up to a cutoff scale Q 

f(p < Q) ~ a" c , /(p > Q) < (Q/p) 4 (2.18) 

so that the high momentum particles (p > Q) do not dominate energy density or anything 
else and can be neglected in the discussion. For c > 0, the soft modes are overpopulated 
compared to the thermal ensemble with the same energy density. We expect thermalization 
to proceed by energy cascading to higher p modes with lower occupancy, until the typical 
occupancy reaches 1 and quantum mechanics cuts off further cascading. We will see that if 
c < 1, both the elastic scattering and inverse collinear splitting are of comparable effective- 
ness in transferring energy to higher momenta. For c > 1, the occupancies are so large, that 
even in the weak coupling limit, the thermalization is driven by non-perturbative physics. 
This will be discussed separately in the next subsection. 
The energy density of the system is 

e~ Jd 3 ppf(p)~ a - C Q\ (2.19) 

and the final temperature is T nna i ~ e* ~ a~ c ' 4 Q, so the equilibration cannot occur faster 

than in time t cq ~ a ~ 2 Tfin a i ~ a~ 2+c ' Q~ l , which is the equilibration time for a small 

deviation from thermal equilibrium in a system at temperature T nna i. On the other hand 

systems with high occupancies have Bose-stimulated interaction rates so we also expect it 

should not take longer than this time scale for the system to thermalize. 

Let us now estimate the time scale it takes for a particle with p ~ Q to appreciably 

change its angle or momentum. For the initial conditions of eq. (2.18), the momentum 

diffusion due to elastic scattering is dominated by modes with p ~ Q and is parametrically 

of order 

Mastic ~ o? / /(p)[l+/(p)] ~ a 2 ~ 2c Q 3 , (2.20) 

Jp 

and the time for large-angle change is 

ilargc angle ~ Q 2 /Q ~ a" 2 "*" 2 ^- 1 • (2.21) 

On the other hand, for particles with momentum p <C Q the time scale is shorter by a factor 
of (p/Q) 2 - Such "soft" particles of low momentum will equilibrate quickly into an infrared 
tail with occupancy f(p <C Q) oc p _1 - the functional form which is in quasi-equilibrium 



10 



in the sense that the "gain" and "loss" terms in the Boltzmann equation cancel. Such 
a soft tail does not dominate either screening, elastic scattering, or the energy density of 
the system, so it plays no role in directing the dynamics of the typical excitations. Such 
f(p) oc 1/p tails will be a common feature in the sections which follow, so we introduce some 
notation to describe them. In our case, at and after the time £i arge angle the particles have 
had time to organize themselves into such a 1/p tail. The distribution is characterized by 
two scales; the maximum momentum p max below which the occupancy behaves as oc 1/p, 
and an "effective temperature" determining the occupancy below p ma , x : 

/(p)~— eow-p). (2.22) 

Of course the Heaviside function is an oversimplification of the transition from oc T*/p 
behavior to rapid falloff; but for parametric estimates it will be sufficient. The parameter 
T* is solved from energy conservation 

e ~ T^Lx ■ (2-23) 

The subsequent equilibration is then controlled by the evolution of the scale p max . Before 
thermalization p max < ^finai < T*, and when p max increases T* decreases. When p max ~ 
Tfinai ~ T*, the typical occupancies at scale p max reach 1, and all scales below Tg na i have 
thermalized. After that, the cascade will continue (without Bose stimulation) to higher 
momentum scales. However, in thermal ensemble the energy density, screening, and elastic 
scattering are dominated by scale T nna i, so that even though the ultraviolet tail of the 
distribution has not completely adjusted to the Maxwell-Boltzmann form, the system can 
be considered thermalized. 

The scale p max can grow either by elastic scatterings or inverse splitting processes. The 
rate for large angle elastic collisions is 

Elastic ~ a 2 T^p mSbK , (2.24) 

y, ^Clastic 2 rri ( ^* \ ( 00 r\ 

1 large angle ~ — ~ « ±* I I , (2.25) 

Pmax VPmax/ 

and is enhanced compared to the equilibrium one ~ c?T . 

The rate for a hard joining process, in which two particles with p ~ p max merge, so that 
the momentum of the resulting particle has changed by order 1, is given by the eq. (2.17) 
with p ~ k ~ (p—k) ~ Pmax so that 

dtf(Pma. x ) ~ r har d merge/(Pmax) ~ ar c l/(p m ax) [l+/(Pmax)], (2.26) 

and inserting the elastic scattering rate, 

2 m I 1* 



Thard merge ~ « % ( J . (2.27) 

The rates for large momentum transfer either from elastic scattering or hard joining are the 
same (up to logarithms) so both the two processes compete, but all parametric estimates 
made using either one will give the same results. 
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Figure 4. Evolution of the momentum scales in a log-log plot during the cascade with initial 
distribution of eq. (2.18) with c = 1. 

During the evolution of p max the total particle number changes, and we have assumed 
that the number changing processes can keep up with the change of p max and hence T*. 
We expect that this is true because number-changing 2 O 1 splitting processes are para- 
metrically as efficient as elastic scattering. But even if they are numerically less efficient, 
particle number is still efficiently changed, because the particle number changing rate is 
large in the infrared. Specifically, the total particle number changing rate is of order 



i"l -H-2" 



dr BH 

dk-g±[l+f(k)}~aT el f(m), 



(2.28) 



which is greater than the hard merging or the large angle elastic scattering rate by a power 

of (Pmax/m) ~ y'pma.x/aT* [19]. 

The time evolution of p max and T* can be determined by self-consistently solving 



Pmax r ^ J (/elastic^) 

together with eq. (2.24) and energy conservation, eq. (2.19) and eq. (2.23): 



Pe 



2-2e 8 1 
~ O 7 Q7tT 



Ti 



-6-c ^4 _3 

~a 7 Qit 7 



(2.29) 



(2.30) 



2c — 2 _-, 

These expressions hold for all times t such that p max > Q, that is, t > a 7 Q~ , until 
equilibration is complete, which occurs when p max ~ T*, at time 



i eq ~ a 2+ iQ x ~ a 2 T final 



(2.31) 



These results are summarized in Fig. 4. 

2.3 Nielsen-Olesen Instability 

In the last subsection we considered the case f(p ~ Q) ~ a~ c with < c < 1. Now we 
briefly comment on what happens in the case c > 1, that is, the case of extremely high 
initial occupancy. The physics in this regime is nonperturbative; large occupation numbers 
of order > a~ 1 can cancel the suppression of interactions coming from the small coupling 
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constant. Specifically, the nonperturbative physics of the Nielsen-Olesen instability [13] 
will dominate in this regime. 

Very large occupancy fields can be viewed as effectively classical. In our case, we 
have classical color-magnetic fields of coherence length £ co h ~ Q^ 1 and of field strength 
B 2 ~ a~ c Q 4 — > B ~ a~ c ' 2 Q 2 . The Lamor radius of a particle propagating in such a 
field is L ~ a~ 1 ' 4 B~ 1 ' 2 ~ cr c " Q , which is shorter than the coherence length of 
the field. Therefore we can neglect the spatial variation of the magnetic field. In an 
intense and uniform magnetic field, relativistic excitations split up into Landau levels with 
energy spacing ~ L _1 . But as pointed out by Nielsen and Olesen [13], the way the spin- 
magnetic interaction splits the Landau levels means that, for spin-1 particles, one spin 
state of the lowest Landau level is actually exponentially unstable, with a growth time 
7 ~ L" 1 ~ a~ ~Q. In a time scale 5 t ~ L^ 1 these unstable modes grow to absorb the 
energy density of the magnetic fields, transferring it to modes of typical wave number 
p ~ L _1 ~ a~^Q 3> Q. So the instability transfers energy towards ultraviolet and works 
towards thermalizing the distribution. Using energy conservation, a~ c Q A ~ fncw{p)p 4 , we 
find the final occupancy of these modes to be f(p) ~ a" 1 . Therefore the final state has 



l-c 



Qncw ~ a 4 Q , / ncw ~ a l . (2.32) 

So the end product of the Nielsen-Olesen instability is a distribution with c = 1. The 
equilibration of this "final" state will then proceed as described in the previous subsection. 

2.4 Landau-Pomeranchuk-Migdal suppression 

We will turn next to dilute systems, c < or f(p) < 1. In this case the screening scale 
m 2 gets smaller, but surprisingly the elastic scattering rate does not; while there are fewer 
scattering "targets," this is compensated by the weaker screening, so the explicit factor of 
f(p) in eq. (2.15) is canceled by the implicit factor in 1/m 2 . But the individual elastic 
scatterings are by smaller and smaller angles. If the scattering is by a sufficiently small an- 
gle, any emitted radiation stays "on top" of the emitter for so long that it can still interfere 
with an emission from the next scattering, which is called the Landau-Pomeranchuk-Migdal 
(LPM) effect [14], after the scientists who understood the effect in the context of QED. 
This effect was first understood correctly in the context of QCD by Baier et al ("BDMPS") 
[15]. These interference effects will prove to play an important role in the case of a dilute 
(low initial occupancy) system, so we will review the physics of the LPM effect before we 
discuss the low-occupancy case. 

Consider the quantum mechanical formation time £f or m of a nearly collinear splitting 
process, defined as the time it takes to separate the wave packet of the emitted daughter 
particle from the emitting mother particle in the transverse direction. Take k to be the 
momentum of the daughter particle and k± its transverse momentum with respect to 
the mother particle. The transverse size of the wave packet is then Ax± ~ l/k±, and 



J The timescale also depends on the size of the initial seed fluctuations in p ~ L _1 modes. But at 
minimum there are vacuum fluctuations, which behave as f(p) ~ 1/2. So the time for the fluctuations to 
grow from / ~ 1 to / ~ a~ 1 is ~ L _1 ln(a _1 ), longer than L~ l by a logarithm. 
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Figure 5. Interference between diagrams in which an emission occurs either before or after multiple 
scattering events. When they occur on the final states, the scattering events can involve either 
external leg. These interference effects are responsible for the LPM effect. 



its transverse velocity v± ~ k±/k, so that the daughter's wave packet overlaps with the 

mother particle until the time 

Ax± 



tfo 



V± 



k_ 



(2.33) 



The daughter particle acquires transverse momentum by random elastic scatterings 
with the other particles in the medium. If tf orm is longer than the mean time between such 
scatterings, we can replace individual scatterings with transverse momentum diffusion and 
estimate that k 2 , ~ <7 e lasticiform> giving the formation time 6 



Cformv'v 



k 



^clastic 



(2.34) 



If the formation time is longer than the time between elastic collisions T^J , the kinetic 
theory treatment of eq. (2.3) is no longer adequate as the diagrams depicted in Fig. 5 
interfere coherently, an effect not taken into account in the kinetic treatment. Physically, 
the emitting particle can no longer resolve individual kicks originating from individual 
collisions, and effectively sees only the net deflection from all collisions during the formation 
time; the particle's trajectory bends smoothly. Then the rate for splitting processes reads 7 



^LPM 



^(k)^at^ r Jk), 



r s piit("0 



min[T s B p H t 



(k),T^(k)}. 



(2.35) 
(2.36) 



Equivalently, the spectrum of split particles changes from Bethe-Heitler log (A;) to LPM 
type k~ l ' 2 spectrum at the momentum scale &lpm where the formation time is of the same 



6 In nonabelian theories, both the mother and daughter particles are subject to elastic scattering as 
the gauge bosons are charged. In an abelian theory the photon is not charged and the formation time 
depends on the momentum of the mother particle instead of that of the daughter. The analogous estimate 
in the abelian case is tf OTnl (p) ~ y/p 2 /kq c \ ast ic, where p,k are the momenta of the mother and daughter 
respectively. 

7 Note that the formation time of an inverse splitting process in the previous subsection (overoccupied 
case) is tform ~ (aT*)~ ■ This is of the same order of magnitude as the time between elastic collisions 



-i 



- el 



(aT*) , so the LPM suppression does not affect the parametric estimates. 
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order of magnitude as the time between elastic scatterings: 

&LPM ^ f m 2 \ 

^clastic \ ^elastic / 

4 

&LPM ~ • (2.37) 

^clastic 

2.5 c < 0: Low occupancy and "Bottom-up" thermalization 

Next we consider the equilibration of an isotropic homogeneous system, whose initial dis- 
tribution is given again by eq. (2.18) with typical occupancies a~ c < 1 (c < 0), so that the 
initial distribution consists of underoccupied hard particles at a typical scale Q > Tfi na i. 
The thermalization proceeds in close analogy to the "bottom-up" scenario described in 
Ref. [12]: First, a population of soft modes with momenta of order ~ m is generated by 
soft splitting, induced by small angle elastic collisions between the hard particles. The 
soft sector subsequently thermalizes in a time Qt ~ a~ 2 ~ c < 3 . By this time, the soft sector 
carries only a small portion of the total energy of the system, but it draws energy from 
the hard modes by hard LPM suppressed splitting. The system becomes fully thermalized 
when the hard sector has lost all its energy to the soft sector at Qt ~ a~ 2 a~ 3c ' 8 . 

2.5.1 Qt < a' 2 ' 

For the underoccupied initial condition, the screening scale, momentum diffusion constant 
due to elastic scattering, and small angle elastic scattering rates read 

Elastic ~ a 2 n h ~ a 2 ~ c Q 3 , (2.38) 

m\ ~ an h /Q ~ a l ~ c Q 2 , (2.39) 

T el ~ ^i£ „ aQ t (2.40) 

m h 

where n^ ~ a~ c Q 3 is the number of initial hard particles. The index h indicates that 

the screening is due to the original hard modes. If the equilibration would proceed via 

elastic scatterings as in the overoccupied case, the estimate for the equilibration time would 

be dictated by the large angle collision time for the hardest modes t cq ~ Q 2 /q c iastic ~ 

a~ 2+c Q~ 1 . However, as we will see, inelastic scatterings provide a faster route to the 

thermal distribution. 

The hard particles radiate soft gluons creating a new population of particles at small 

momentum. The production rate of gluons with soft momentum k by hard particles with 

momentum p is Bose stimulated. The production rate of soft particles of momentum k, 

from eq. (2.17), is 



r pro d(A0 ~ / r sp i it (A;)/(p)[i+/ s (fc)][i+/(p)] ~ / r split (£0/(p)[i+/ s (fc)] , (2.41) 
jp jp 

where f s is the distribution of the new soft modes. But the stimulation factor [l+f s (k)] is 
cancelled by the inverse process where the soft gluon merges onto a hard particle 

Tabsorb ~ f r split (k)f s (k)f( P )[l+f(p)] - [T split (k)f s (k)f( P ) , (2.42) 

J p Jp 
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so that the total production rate of soft particles becomes 

(r pro d - Tabsorb) ~ / T sp n t (k)f(p) . (2.43) 

Jp 
The spectrum of the emitted particles is then 

/ s (p)~n h r split (p)t/p 3 , (2.44) 

where p~ s arises from the momentum space element. Using eq. (2.16) and eq. (2.35), 

f s (p) ~ an h r ci t/p 3 , for p < k LPM , (2.45) 

f s {p) ~ an h s/4^~t/p 7/2 , for p > k LPM , (2.46) 

with (see eq. (2.37), eq. (2.38) and eq. (2.39)) 

&lpm ~ oT c Q. (2.47) 

As long as c > — 1 so that /clpm > TUh , the production of the softest modes is described by 
the Bethe-Heitler rate. 8 However, interactions can change the actual spectrum from this 
production rate. In particular, an f(p) ~ p~ 3 spectrum rises more steeply than a thermal 
spectrum. Below some scale p max the spectrum will collapse into a thermal-like form, 

fsip) ~ T m /p, for p < p max (2.48) 

in close analogy to the overoccupied case discussed in subsection 2.2. The soft sector's 
contribution to the screening and elastic scattering are dominated by the scale p ma ,x, as 
they were in the overoccupied case. The major difference is, however, that the growth of 
the soft sector is driven by soft splittings from the hard particles. In addition to interacting 
among themselves, the soft particles can scatter with the hard particles. 

Let us now determine what physics dominates the evolution of p m ax- There are four 
mechanisms which compete in redistributing the soft sector: 

• Elastic scattering with hard particles, described by ^clastic : The scale which has had 
enough time to experience order 1 momentum changes is 

Pmax ^ ~ (^elastic*) 172 ~ OC^^QiQt) 1 ' 2 . (2.49) 

All scales below p^l^ scat have had time to redistribute, via elastic scattering, into 
an f(p) oc 1/p thermal-like tail. 

• Elastic scattering among soft particles: That the integral in eq. (2.13) diverges for 
oc 1/p 3 already shows the inconsistency of such a soft distribution. While the total 
number of soft particles is initially smaller than that of the hard particles, their 



For extremely dilute initial conditions c < — 1, the production of even the softest modes is LPM 
suppressed. In this case the first soft gluons emerge at the time t ~ tf orm (mh) ~ a~ 3 ' 4+c ' 4 Q~ 1 , after which 
the evolution will follow the normal c > — 1 case discussed in the next sub-subsection. 
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interaction rates are Bose stimulated, so it could be possible that elastic scattering 
is dominated by the soft sector. 

Assuming that the evolution of p max is dominated by the elastic scatterings between 
soft particles, we can estimate the momentum diffusion constant arising from the 
soft sector (^elastic ~ o> 2 T*p^* c 8cat . This momentum diffusion coefficient re-arranges 
the soft particles out to p max ~ (<Zs,eiastic^) 5 - The energy density of the soft sector is 
dominated by the particles emitted at scale p max , so the energy density of the soft 
sector is e s ~ n/ i r sp i it (p^ x scat )£ ~ T*p;^ ax . Combining and solving self-consistently, 
we get 

p softscat _ a 6/5-2c/5Q(Q t) 3/5 _ ^ ^ 

Merging at the p max scale: The particle number in the overoccupied regime can 
decrease by merging pairs of particles to create fewer particles at higher momentum, 
analogously to the hard merging in the overoccupied case. All modes up to a limiting 
scale Pmlx mcrge have had time to undergo repeated merging to push the extra particle 

i ■. ,i i hard merge 

number to the scale p m ax • 

To estimate the scale pmax merge we find the hardest modes which have had time to 
have an order 1 probability of merging 

rmerge^S 1116 ^)*-! (2.51) 

with the Bose stimulated rate 

rVrgeCp) ~ r split (p)[i+/(p)] , (2.52) 

so that 

p hard merge _ Q 4/3-c/3g(Q t )2/3 _ (2.53) 

Saturation: For a given soft momentum i<Q, emissions will cause f s (k) to grow 
until it reaches a size where absorption (joining) processes shut off the production. 
Let us find the occupancy where this occurs. The production and absorption rates 
per hard particle f(p) are 

f{ P ) (r prod - r absorb ) ~r,^ Iit (fc)(/(p)[i+/(p-fc)][i+/ a (fc)] - /(p-fc)/.(fc)[i+/(p)] 

~r split (fe) (kf a (k)£ + f(p)[l+f(p)]\ . (2.54) 

This rate, integrated over p, goes to zero when 

m ~ J,-**/* ■ ( 55) 

Since the phase space opens up as p 2 dp, it must be that df(p)/dp is on average 
negative. So we can estimate df /dp ~ —f(p)/p- Using this estimate and using that 
p ~ Q is the dominant scale, we find 

/ 8 (fc)~|[l+/(p)]~|(l + a- c ). (2.56) 
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(The 1 + a c term is included so the formula can also be applied for overoccupancy.) 

Whenever the occupancy of soft modes reaches this value, soft modes are absorbed as 
quickly as they are emitted. The production mechanism shuts off and the occupancy 
grows no further. If this shutoff is responsible for deciding where the momentum tail 
goes from k~ 3 to k~ l behavior, then the scale p max is determined by the condition that 
the emitted spectrum, / s (p max ) from eq. (2.44), has / s (p max ) ~ Q/Pmax- Inserting 
eq. (2.45), we find 

^saturation _ a l-c/2 Q ^l/2^ (2 ^ 

which turns out to be the same as the scale set by elastic scatterings with hard 
particles. 

At early times the hierarchy of these scales is 

hard scat saturation ._ hard merge ._ ^soft scat /<-> ro\ 

"max rsJ "max -^ "max -^ "max ^Z.uoJ 

so that the evolution of p max is set by elastic scatterings with the hard particles (as well as 
saturation of soft particle production): 



Pmax ~ PS "** ~ a^QiQt)* . (2.59) 

The occupancies in the infrared tail can be solved again by estimating how much energy 
density the soft sector has had time to eat. The energy density in the infrared tail is 
dominated by the scale p max so that 



E s ~ r„j4ax ~ Pmax^sfrmax) , (2.60) 

BH 



4 

n s(Pmax) ~ r sp i it {PmaxjtlT'h i 

split 



r* ~ if* (p max )tn h /p 2 max ~ Q, (2.61) 



showing that T* is a constant scale. 

There are a number of assumptions we have made whose breaking may lead to a 
revision of the time evolution of p max and T*. First, we assumed that the hard sector 
dominates elastic scattering, screening and energy density. Of these, screening is most 
sensitive to soft modes. The contribution of the soft sector to screening is dominated by 
the scale p max , 

9 f fs(p) ^ nTn __ _ 2-c/2 n 2/ n ,Nl/2 



m ' I i^L „ a r,p max ~ a^Q^QW, (2.62) 

J P P 

which becomes comparable to hard screening at the time (compare eq. (2.62) to eq. (2.39)) 

Qt ~ a' 2 ' . (2.63) 

Second, we also made the assumption of the hierarchy given in eq. (2.58), which also breaks 
down at the same time Qt ~ a~ 2_c , when all four scales reach 

hard scat saturation hard merge soft scat , ( 2 64) 

fraax rmax max rmax ^LfM- V-" - "^/ 
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2.5.2 a~ 2 - c <Qt< Q" 2+c / 3 

In the second stage, the scale p max continues its growth, but as it enters the LPM suppressed 
regime, the increase of energy to the soft sector is no longer enough to keep T* a constant 
scale. At the same time, the soft sector grows large enough to dominate screening, so that 
the rate of elastic scattering reduces and the new LPM scale follows p max . 

Assuming that the evolution of p max is still set by the elastic scattering off the hard 
particles, we get the same estimate as before, 

Pmax ~ (^clastic*) 172 ~ a 1_c/3 Q(Qt)5 . (2.65) 

Then, from the energy density that the infrared tail has had time to eat we get 

T* ~ T};™(p max )tn h /p 2 max ~ a^-^iQty^Q. (2.66) 

The contribution of the soft screening and elastic scattering from soft particles are 
then 

m 2 „ a 3/2-3c/4 Q( Q f)1 /4 j (2 67) 

^.elastic ~ a 2 ~ c Q 3 , (2.68) 

so that the soft and hard sectors give comparable contributions to elastic scattering but the 
soft sector dominates screening. Hence, the assumption that p max can be estimated from 
hard elastic scattering is selfconsistent. Furthermore, as the dynamics of the soft sector 
is set by the soft sector itself, we know already from the overoccupied case that the hard 
merging scale will equal p max . 

The growth of m s decreases the elastic scattering rate, and so the LPM scale starts to 
grow. The new LPM scale is then 

&LPM ~ "4/<?clastic ~ Pmax, (2.69) 

so that using either the BH rate or the LPM rate will give the correct estimate for the 
energy density of the infrared tail. 

When p max ~ T*, the occupancies at the scale p max become of order 1 and the soft 
sector thermalizes. This happens at Qt ~ a _2+c ' 3 , and at later times, the estimates again 
need to be revised. 

2.5.3 a~ 2+c / 3 <Qt< a~ 2+3c / 8 

In the final stage of the equilibration the soft particles dominate screening, particle number 
and q but the hard sector still carries most of the energy. The soft sector is thermalized 
and is described by a single scale, its temperature T, with /clpm ~ T, n s ~ T 3 , e s ~ T 4 , 
and 

^elastic ~ a 2 T 3 . (2.70) 

The hard sector loses its energy to the soft sector by splitting off daughters, which split 
further, cascading down in momentum to join the thermal bath. There is a subtlety, first 
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identified by Baier et al [12] . As discussed, the spectrum of excitations emitted by a hard 
p ~ Q particle scales as k~ 1 ' 2 dk/k. In terms of energy lost, the spectrum is dominated by 
the largest k; in terms of number of excitations, it is dominated by small k. But in a given 
amount of time t, a sufficiently high-energy daughter will not actually give its energy to 
the thermal bath, but will remain intact. In terms of energy gain by the thermal bath, the 
most important emitted daughters are therefore the highest-energy daughters which can 
break down completely and yield their energy to the thermal bath within time t. These 
are precisely those excitations which can hard-split in less than time t. The rate for a 
particle of momentum k to hard-split is of order the same as the rate for a much harder 
p ~ Q particle to emit a particle of momentum k. Therefore the most important daughter 
momentum scale k is that, such that each hard particle emits 0(1) such daughters in time 
t. Call this momentum scale fe s piitj parametrically (see eq. (2.34), eq. (2.35) and eq. (2.70)) 

r^f (k sp m)t ~ 1, fc split ~ a 4 Th 2 . (2.71) 

The particles with momentum fc sp iit then undergo multiple splittings, in a time scale which 

is much shorter than the formation time of the initial emission, depositing all their energy 

to the soft thermal bath. The energy density and the temperature of the soft sector are 

then 

i 
e s ~ n h k spUt , T ~ el , (2.72) 

and solving self consistently we get for the temperature and for the splitting scale 

T ~ a 4 n h t 2 ~ a A - c Qh 2 , (2.73) 

fcsplit ~ a 16 ~ 3c Q(Qt) s . (2.74) 

By the time the splitting scale reaches Q, all of the hard modes have had time to 
undergo hard splitting and join the thermal bath, and the system thermalizes when 

Qt~a- 2+ ' ic '\ T final ~a" c / 4 Q. (2.75) 

The final equilibration time should be compared to the characteristic time scale t ~ 
a~ 2 T ~ a~ 2+c > 4 Q for the thermal bath to self-equilibrate. The equilibration time is 
longer but only by a factor of \JQJT ~ a c < 8 ~ /initial- ^ n particular, equilibration is 
much faster than the naive large-angle scattering time for the original hard excitations, 
t eq ~ a~ 2+c Q~ 1 . Physically, the factor \JQJT arises because the time it takes a particle 
of momentum Q ^> T to radiate away its energy in a thermal bath of temperature T is 
t ~ (a 2 T)~ 1 ^Q/T. So the final equilibration time is just the time for the hard excitations 
to break up in a thermal bath at the final temperature. 

3 Weakly anisotropic systems 

In some applications we expect a nonabelian plasma to form, in which, rather than being 
evenly distributed in direction, particle occupancy is larger in some directions than in 
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Figure 6. Evolution of the momentum scales in a log-log plot with initial distribution of eq. (2.18) 
with c = — 1 . 

others. In some cases the excess occupation may be modest. For instance, a nonabelian 
plasma which is experiencing shear flow develops an anisotropic particle distribution. This 
has been discussed by Asakawa et al [20], who argue that in a range of shear flow strengths, 
the plasma remains nearly isotropic but only because of the physics of plasma instabilities. 
Weakly anisotropic plasmas should also be rather generic in the presence of phase interfaces 
during early Universe phase transitions, such as the electroweak phase transition (if it 
exists). 

After an overview of the physics of plasma instabilities, we will first treat this case 
of weakly anisotropic plasmas, and then consider plasmas with a very high degree of 
anisotropy in the next section. Since anisotropic plasmas are rather complex and because 
the "parameter space" describing them is larger than for the isotropic case, we will not 
attempt to follow the dynamics through every stage to equilibrium. Instead we will iden- 
tify the relevant dynamics initially, and then determine the time scale before the dynamics 
change appreciably and what the new relevant dynamics becomes. 

3.1 Plasma instabilities 

So far we have only considered systems which are statistically locally isotropic. A new com- 
plication enters when we consider systems which are locally anisotropic. Namely, certain 
long-wavelength magnetic gauge field excitations are then generically unstable to expo- 
nential growth, in a process called the Weibel instability (see [6] for the electromagnetic 
case and [7-10] for an overview of the nonabelian case and [21-23] for some nonabelian 
numerical studies). At weak coupling, in many circumstances these fields are expected to 
grow large enough that they come to dominate much of the dynamics of the system. In 
particular, the "hard" excitations which dominate the energy density of the system may 
predominantly change direction by deflection in these magnetic fields. These deflections 
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may then dominate the splitting or joining processes these excitations undergo. The key 
physics involved in plasma instabilities is the physics of screening. So far we have only 
needed a very crude description of plasma screening, but plasma instabilities involve finer 
details of the screening process. Therefore we begin with a more careful discussion of the 
physics of screening. 

There are already nice, physical and quantitative discussions of plasma instabilities in 
the literature, see for instance section II of Ref. [9] and Refs. [8, 10]. Nevertheless we will 
give a quick qualitative overview to capture certain salient points for our discussion. We 
start by explaining how to think about the screening scale m we previously introduced. 
Consider the case where many particle- like excitations 9 all move along the ±x axis, in the 
presence of a magnetic field with wave-vector k in the y direction and B field in the z 
direction, as illustrated in Fig. 7. The force on a particle is 

F = q a v-B a . (3.1) 

If the particle trajectories are undeflected at time t = 0, then by time t their velocities and 
positions have changed, to linear order in B, by 



q a B a t . . , q a B a t 2 

p 2p 



%(*) = — ~ > Mi) = ^TT^ • (3-2) 



This particle motion concentrates particles in some regions and depletes them in others. 
The concentration fraction is kSy{t), and the net current associated with this concentration 

is 



J a ~ J d*p f(p)q a k5y(t) ~ J d*p f(p)^L k B b t 2 . (3.3) 

Here J d 3 pf(p) is just counting the density of particles. We recognize the combination 

rf f , n9V x 2m fd 3 pf(p) 2 /„ „N 

dp f(p) = Sab 9 T r / 5 ab m (3.4) 

P J P 

so the induced current is 

J a ~ kB a m 2 t 2 . (3.5) 

This current enters the equations of motion for the magnetic field, where it competes with 
terms of form VxB~ kB. Therefore the current term becomes important in the evolution 
of a magnetic field whenever m 2 t 2 ~ 1. So the correct physical interpretation of the thermal 
mass m is that it is the inverse of the time scale for particle deflections to add up to enough 
that they influence the gauge field dynamics. That is, if particles remain coherently in a 
region of electric or magnetic field of one sign for time scales of order \/m or longer, then 



9 By particle-like excitations we mean, excitations which can be described in terms of wave packets which 
are of smaller extent than the coherence length of the classical fields in which they move; and yet are built 
of a narrow enough distribution of momenta that they move nearly non-dispersively. If the typical particle 
momentum is p, the range of momenta particles carry (and therefore the range of momenta out of which 
wave packets are built) is Ap, and the characteristic wave number of the background fields they move in is 
k, this requires p S> Ap S> k. We will encounter a case where these criteria cannot be met in subsection 
4.8, but in most of our discussion the scale ordering will be true parametrically. 
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Figure 7. Illustration of the Weibel instability. Particles flowing in the z direction are deflected by 
a magnetic field with both polarization and wave- vector at right angles to the particle flux. Left: the 
deflection of positive charged particles (in blue) focuses them into some regions and dcfocuscs them 
from others. Right: opposite charged particles (red) are focused in different regions but generate 
the same-sign net currents. Both are focused such that their currents amplify the original magnetic 
field. 

the induced currents are large enough to significantly influence the field evolution. In the 
case illustrated, since the particles enhance the magnetic field which deflects them, they 
lead to an exponential growth in the strength of the magnetic field B(t) ~ B(t = 0) exp(7i), 
with a growth rate 7 ~ m. 

For the case considered in Fig. 7 the currents from the particle excitations strengthen 
the magnetic field, which destabilizes the plasma towards magnetic field growth. Particles 
flying in certain other directions contribute currents which oppose the magnetic field. For 
an isotropic distribution of particle momenta, a magnetic field will not induce any current. 
That is because a magnetic field induces an overall rotation of the distribution of particles; 
but if the particle distribution is isotropic then a rotation does not change it. Therefore 
an isotropic distribution is neutral in the sense that it neither stabilizes nor destabilizes 
magnetic fields. For an anisotropic distribution of particles, we can instead ask about 
the average over directions and polarizations of magnetic fields. This is equivalent to 
averaging over the particles' momentum distribution. Hence, for anisotropic systems, the 
plasma's impact on magnetic fields, averaged over directions and polarizations of magnetic 
fields, is neutral. That means that there will be some directions and polarizations where 
the particles stabilize the magnetic fields (generate currents which oppose the magnetic 
field), while for other directions and polarizations they will destabilize the magnetic fields 
(generate currents which enhance the magnetic field). Hence, in an anisotropic plasma, 
magnetic instabilities will always be present in some field directions and polarizations [9]. 

In practice we will be interested in magnetic fields which change with time, which 
always means there are also electric fields present. Electric fields are stabilized by an 
isotropic distribution of particles, so time-changing fields tend to be more stable than 
static magnetic ones. Also, if the electromagnetic fields vary on too short of a length scale, 
particles will not have trajectories which stay in a coherent region of field strength for the 
requisite t > 1/m needed for currents to dominate the Yang-Mills equations of motion. 
To make somewhat more precise statements about which modes are unstable under which 
circumstances, we will sketch a diagrammatic evaluation. 
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At the linearized level, an instability will occur if the retarded propagator has a pole for 
some wave number k and pure imaginary frequency u = ij, corresponding to exponential 
growth. The transverse part of the inverse retarded propagator is 



G#(ft,7) 



+ l 2 )ik 



K% Kj J 



IIij(fc,7) 



(3.6) 



If the self-energy loop momentum p is large compared to the external momentum k, then 
the self-energy behavior becomes spin independent and we can treat the simplest case, 
which is a scalar in the loop. For this case, 



lF"(fc) ~ g 2 / d 4 p 
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(3.7) 



(p — k) 2 

Assuming p ■ k ^> k 2 , we can expand in 2p ■ k 3> k 2 . Together with the spin-independent 
form for the self-energy, this constitutes the hard-loop approximation for the self-energy 
In this approximation, the self-energy is 
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(3.8) 



Positive-sign contributions inside the square brackets contribute positively to eq. (3.6), 
stabilizing the field. Negative contributions destabilize the field, and if they are large 
enough - if the RHS of eq. (3.8) equals or is larger than 2(fc 2 + 7 2 ) - then they allow an 
exponential instability. We see that those particles in the angular range v ■ k < 7 are the 
ones contributing to the instability; all others stabilize. 

The above discussion is based on linearizing in the size of the magnetic field, and 
possibly in other interactions. It breaks down when higher loop contributions to the self- 
energy become as large as the one-loop self-energy discussed above. Most of the time, the 
dynamics controlling the equilibration of the full system take place on a time scale long 
compared to I/7 the growth rate of the instabilities. In this case, the instabilities will grow 
until some nonlinear physics limits their size. After this, the relevant magnetic fields will 
remain large, evolving but in a "quasi-stationary" fashion so that statistical properties of 
their distribution evolve only on the time scale of the system's approach to equilibrium. The 
saturation of instabilities and quasi-stationary evolution have been verified in numerical 
studies [11, 24, 25]. 

3.2 Instabilities for weak anisotropy 

Consider a system along the lines of what we discussed in the previous sections; that is, 
a distribution of "particles" with a typical momentum scale p ~ Q, sufficiently fast falloff 
in f(p) above this scale, and typical occupancies f(p) ~ a~ c for p ~ Q. We will consider 
both the possibility c > (overoccupancy) and c < (underoccupancy). But now we will 
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consider the case where f(p) is also a function of the direction p. In this section we will 
take this angular dependence to be weak, 

fip) = f(p) x (1 + eF(p)) (3.9) 

with F(p) some order-1 combination of spherical harmonics of number 2 and higher, but 
dominated by low numbers, 10 and e characterizing how small the non-spherical component 
of f(p) is. We also define d as e = a (note that d is negative) to make it easier to mix 
powers of a and of e. 

Inevitably the distribution we consider here is not the most general we could consider. 
While there is no need to allow Y\ m (p) content in F(p) we could have considered cases 
where f(p) does not factorize into a radial and an angular function, that is, where the 
typical particle energy is larger in some directions than in others. 

First we compute the screening scale and elastic scattering q, 



m 2 ~ aQ 2 f(Q) ~ a^Q 2 , q elastic ~ a 2 Q 3 /[l+/] ~ { 2 _,^ ' (3.10) 



To determine whether plasma instabilities play any role in this system, we need to determine 
what gauge field modes become plasma-unstable, how large the amplitudes of the associated 
magnetic fields grow, and how much deflection these magnetic fields cause in the hard 
particles. That is, we need to compute q\ ns t arising from soft-unstable magnetic fields and 
compare it to (/elastic to determine whether the instabilities play an important role. 

First let us see how many modes are unstable and how fast they grow. As discussed in 
the last section, a mode is unstable and grows if the inverse retarded propagator, eq. (3.6), 
vanishes for an imaginary frequency u = ij. —Hij gets more positive as 7 is increased, 
so the range of k which are unstable is determined by considering 7 = + . The isotropic 
part of the integral vanishes for 7 = + and gives ivy/k for small 7. The anisotropic part 
averages to zero over k directions, but is positive in some directions and negative in other 
directions, e.g. in unstable directions 

„ m 2 /7T7 \ , 

- n «~-(T- £ )- (3 ' u) 

Applying this we find instabilities in about half of directions, with maximum k of 

ki st ~em 2 . (3.12) 

A simpler way to get this scale would be to compute m 2 just from the anisotropic piece 
of the particle distribution, which gives fn 2 nst ~ em 2 . The unstable modes are those with 
kf nst ~ fifnsf So long as e < 1 and c < 1, the scale h 1DS t <C Q and so the "hard-loop" 
approximation used to establish these estimates is self-consistent. 

The maximal growth rate is roughly 7 such that the m 2 Ti^/2k term cancels the negative 
term, giving 

3 

7 ~ e^'inst ~ earn. (3.13) 



1 That is, the Yoo(p) component defines f(p). Any Y\ m (p) component can be removed by a choice of 
rest frame. 
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Figure 8. A cartoon illustrating how color rotation could inhibit plasma instabilities. A bluc- 
anti-red magnetic field deflects blue and red color charges. However where color rotation occurs, 
the induced current will be in another color direction (such as green). This does not amplify the 
inducing magnetic field because it is in the wrong color direction; hence plasma instabilities are 
inhibited. 

This is not the answer we would get by just considering the contribution of the unstable 
modes. That is because a time-varying gauge field always has an electric component, and 
an isotropic distribution of particles exerts a stabilizing effect on electric fields. So while 
the isotropic part of the particle distribution does not change which modes are unstable, 
it can slow down the growth of the unstable modes. (This is an example of the general 
phenomenon that scales k <C m have slow dynamics [26]). 

This analysis of which modes are unstable is perturbative. The mode amplitudes grow 
until the perturbative treatment breaks down, which will occur when some nonabelian term 
becomes important somewhere in the above calculation. The evaluation of the hard loop 
assumes particles can propagate a distance \/h ins t without color rotating. Color rotation 
destroys the growth of plasma instabilities, as illustrated in figure 8: when particles color- 
rotate, the current they induce is in a different color direction than the magnetic field which 
generated the current. Hence the current does not enhance the generating magnetic field. 
Color rotation is important when a Wilson line U = Pexp ig J dl ■ A of length / ~ fcj~ st 
generically contains an order- 1 phase factor. A gauge- invariant criterion, for such phase 
factors to be large, is that a Wilson loop of size l/Mnst on a side have an order-1 phase factor, 
see Fig. 9. The phase factor of such a Wilson loop is g~ 1 B/kf nst . So color randomization 
is important if B 2 ~ a~ l kf nst . This translates to an occupancy for the unstable modes of 
0(a~ 1 ). In other words, the growth of plasma instabilities is shut off when plasma-unstable 
modes have nonperturbatively large occupancies. We would arrive at the same conclusions 
by investigating where mutual interactions between such modes significantly influence their 
evolution. We would also reach the same conclusions by asking where two-loop self-energy 
corrections grow to be as large as the one-loop self-energy correction we studied. 
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Figure 9. Setup for the Weibel instability as in Fig. 7, but with a Wilson loop indicated. When 
the phase around the Wilson loop is order- 1, then in any gauge one or both of the indicated particle 
trajectories suffers a large color rotation. Choosing the Wilson loop to be of length ~ 1/m and 
width ~ 1/fc, an order-1 phase in the Wilson loop is sufficient to ensure that the type of color 
rotation illustrated in Fig. 8 will occur over the length scale relevant for the development of the 
instability. 

3.3 Momentum diffusion from the instability 

Next we compute the mean squared rate at which these unstable modes deflect the exci- 
tations of momentum ~ Q. Such an excitation feels a time varying force of order F ~ gB. 
This force remains coherently in the same direction for the characteristic coherence length 
of the magnetic field, £ co h ~ k^. (This is also the typical time scale it takes the particle 
to rotate in color, which would also randomize the direction of the force.) Provided c < 1 
and e < 1, the momentum accumulated in one coherence length, 



Ap ~ t coh gB 



k 2 

1,-1 „ " , inst _ 
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Inst' 



1 l-e 

Inst ~ e 2 a 2 Q 



(3.14) 



is small compared to the typical particle momentum Q. Therefore the change in particle 
momentum will accumulate as many small incoherent "kicks" and can be well characterized 
by a momentum diffusion coefficient qi ns t, 



<?inst ~ F 2 t C oh ~ aB 2 /k h 



(3.15) 



We already saw that the magnetic field saturates at an amplitude B 2 ~ a 1 kf nst . Therefore 



Qinst ~ olB i co h 
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(3.16) 



The instability is important if this q dominates the scattering from elastic processes, but 
we expect it to be irrelevant if it is smaller than the scattering already occurring due to 
normal elastic processes. So the domain where the instability is important is 



(Zinst > ^elastic 5 



(3.17) 
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leading to 



{ l + c 
i-c (3.18) 

e > a 3 c > 0. 



When the instability is important, it can change two things. One thing it can do is 
to enhance angle change, erasing e. The other thing it can do is for this angle change to 
enhance splitting/joining processes. This can cause c and e to evolve, and in the case of low 
occupancy it can create a "cloud" of soft excitations which can dominate the anisotropy of 
screening even if they do not replace the hard particles numerically. It cannot lead to an 
evolution of c through elastic scattering because soft unstable fields are primarily magnetic, 
and only electric fields change particle energy. To estimate the electric field subdominance, 
note that E ~ dtA while B ~ d x A, so E ~ (7/fci ns t)-B ~ eB. Therefore the diffusion of 
particle energy ~ E 2 t co h is slow compared to the diffusion of particle direction. The time 
scale for particle energy to change is correspondingly longer than the time scale for e to be 
erased, ending the instability. 

The rate for e to change appreciably due to scattering is the rate for large angle change, 
which is 

tr 1 , ~ — ~ elat (1 " c) Q (3 19) 

''largeanglc q2 ** \<j.-luj 

This time scale is larger (rate smaller) than the growth rate of the unstable modes, 7 ~ 

3 1 — e 

eaa 2 Q, so the treatment is self-consistent in that the unstable modes have time to grow 
before they significantly modify e. 

To find the rate at which splitting/joining processes occur, we need to determine the 
formation time of emitted particles. Moving through the soft-unstable fields, an emitted 
particle picks up a kick of kf nst every length l/Mnstj which is the effective inter-scattering 
distance. The formation time for radiation of a particle of energy k is therefore 

tfarm(fc) ~ p- ~ -. =* tta(i) ~ khr n l ~ fc*e~i a"! ^Q"! (3.20) 

which is always longer than l//c; ns t. Therefore split /join processes are always in the LPM 
regime. The rate of these processes is a factor of a[l+/(A;)] over the formation time; 



r sp iit/join(*0 ~ <x[l+f(k)WoL ~ k-2eia-* c Q2[l+f(k)] . (3.21) 



To determine whether this rate is physically important, we have to consider the physics of 
the particles generated. 

3.4 Joining and e suppression for overoccupancy 

For the case c > (f(Q) > 1) we need not worry about splitting processes; we saw in 
section 2.2 that f(p < Q) scales at most as 1/p, which is too soft to dominate screening. 
We also see, combining eq. (3.21), eq. (3.19) and eq. (3.18), that r sp i it /j oin £i argc ang i e ^C 1. 
Therefore the particles randomize direction, and e is degraded, before particle number can 
change appreciably. 

Nevertheless it is possible that e is significantly affected by joining processes. The rate 
of this process is r split /j oin (Q) ~ eiai^~ c 'Q (eq. (3.21) with k ~ Q and / ~ a~ c ). But what 
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this expression disguises is that this rate is anisotropic in direction, at the order one level. 
The reason is that q arising from plasma instabilities is in general strongly anisotropic, and 
peaked for those particles which are overoccupied. So what can happen is the following. 
Only an 0(e) fraction of particles are anisotropic; but they generate instabilities, which 
cause an order-one anisotropic q. This makes the rate of joining processes faster, by an 
order-one amount, in the directions which have an excess of particles. Joining reduces the 
number of particles moving in those directions, which is precisely a reduction in e. 

The anisotropy is erased by the joining when an e fraction of particles in the overoc- 
cupied direction have had time the join, creating an e change in f(p). This takes place at 
the time 

*e~ £r it/joln~ e4 « = ^ ! 0" 1 - ( 3 - 22 ) 

This must be compared to the rate at which e is degraded by angular deflection; the 
anisotropy is also erased in a time it takes to redistribute the directions of particles 
^large angle- Whichever of these processes is faster dominates, so that the dominant mecha- 
nism bringing down e is particle joining if 

"e "^ '■large angle 

e < o^ (1 - c) . (3.23) 

Otherwise anisotropy is mostly removed by angle change due to q. The time scale for 
isotropization is 

/ t e ~aa^Q~\ e<al( 1 - c ), 

Change ~ < 3 3 (1 _ c) 1 i (1 ) (3.24) 

[ilargeangle~e 2 a 2i >Q l , € > Ot^ 1 c > . 

In either case, after time £ c hange the plasma has become nearly isotropic, and equilibration 
proceeds as described in subsection 2.2. Note that, in every case, Change is short compared 
to oT 2 TI~ ,, the full equilibration time found in subsection 2.2. 

3.5 Soft emissions and e for underoccupancy 

Next consider the case c < 0, so f(Q) < 1 (underoccupancy). In this case the rate of hard 
particles to split off softer particles is larger than the rate for them to join. Soft particles are 
more efficient at screening than hard particles. The splitting rate is also highly anisotropic, 
because q is. Therefore, more soft particles get generated in the directions where there is 
already a particle number excess. This tends to increase the anisotropy of screening, 
potentially enhancing the strength of plasma instabilities. In fact we will see that this 
happens over most, but not all, of the parameter range of interest. 

The soft particle production is in the LPM regime and, analogously to the isotropic 
case, the spectrum of the soft particles is given by eq. (2.46), eq. (2.48) (replacing ^elastic 
with the dominant (ftnst)) 

f s (p) ~ an h ^q- mst t/p 7/2 , for p > p max , (3.25) 

f s {p) ~ T*/p, for p < p max . (3.26) 
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In this case the evolution of p max is dominated by saturation, eq. (2.56), so that f s (p) is 
limited from above by f s (p) ^ Q/p and T* ~ Q. We find p max by equating eq. (3.25), 
eq. (3.26) and using eq. (3.16): 

T* ~ Q, (3.27) 

Pmax ~ ^ t {an h /Qf'h^ ~ e^a^-c)Q { Qtf/?> . (3.28) 

As the gi ns t is 0(1) anisotropic the soft particles are created with an 0(1) anisotropy. 
The modes are then driven towards isotropy by momentum diffusion dominated by g; ns t- 
The anisotropy of the soft modes is described by e s (k), and how anisotropic they are 
depends on the momenta of the soft particles and how large t is. The momentum scale 
which is just becoming isotropic is 

fciso ~ y/E*t ~ ef a^-^QiQt)^ . (3.29) 

At early times p max is larger than k[ so , while for large t the larger of the two scales is h lso . 
In particular, the time when k- lso catches up with p max is 

tiso(Pmax) ~ -^ ~ e-Sa-s^-^Q- 1 . (3.30) 

^inst 

For A; > fci so , the soft particles have the order 1 anisotropy they where created with, 
e s (k > fei S o) ~ 1. Below fcj so , the particles that were created less than ti SO (k) ~ k 2 /q ms t 
ago are still anisotropic; particles created longer ago have had time to undergo large angle 
change and isotropize. The fraction of anisotropic particles, and hence e s , at the scale 
k < k[ so is then 

e s (k) ~ t iso (k)/t ~ A; 2 /4o for A; < A; iso , (3.31) 

e a (fc) - 1 for fc> k iso . (3.32) 

The new anisotropic population of soft particles can subsequently cause new plasma 
instabilities which contribute to gi ns t> an d may eventually dominate it. The contribution 
to q s i ns t from the new soft particles at the scale k is 



Qs, inst 



(*) ~ e 3 J 2 (k)mUk), (3.33) 



where m s {k) is the screening from particles at scale k. The scale p max always dominates 
soft-particle contributions to screening: for k < p max , the distribution only grows as 1/k 
and screening depends on f kdkf(k). On the other hand for k > p max the soft distribution 

7 5 

falls as k~2<Jk, so that screening falls as J k~^dk, and gets its dominant contribution again 
from the scale p max . Then, if p max > h lso so that the particles which dominate screening 
are order-1 anisotropic, then also q s ,mst is dominated by the particles at the scale p max . In 
this case, the new q s ,mst is 

g s ,inst ~ m'(?w) ~ o?l 2 V ^Tl 12 ~ e^a^( 17 " 7c )g 3 (Qi) 3 / 5 , for p max > fc iso . (3.34) 
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On the other hand, if ki so > p max , then the scale p max may dominate screening but it will 
not dominate soft contributions to q as e s (/c) oc k 2 decreases too fast with decreasing k. 
So the contribution to q scales as f lso k~^k 2 dk ~ f lso k~^dk which is dominated by k{ so , 
and 

^,inst ~ m\ (fcto) ~ [aklj s (k iso )] 3/2 ~ a 3 nl /2 q7 n H 8 tV8 (3.35) 

__9_ 39- 15c o 3 

~ e iea i6 Q 3 {Qt)8, for p max < /c iso . (3.36) 

As the soft sector grows, the new gainst increases and it may start to dominate over 
the gi ns t from the hard particles; we define i c hange as the time when 

Qs,mst{tchzmge) ~ <Zinst • (3.37) 

We then have three possibilities: 

• The soft gainst starts to dominate q at early times when particles at p max are still 
fully anisotropic, p m ax(*change) > fti SO (£change)- I n this case the soft sector takes over 
q inst at time 

Change ~e 4 « 4 Q ■ (3.38) 

This is consistent with the assumption for (see eq. (3.30), eq. (3.38)) 

l + c 

^change < iiso(Pmax) => € < a~ . (3.39) 

This is the Region 2 in Fig. 10. 

• The new plasma instabilities start to dominate after the scale p max has isotropized, 
but before time £i ar ge angle which is when k[ so ~ Q and momentum diffusion has 
isotropized all modes. In this case the soft sector begins to dominate g; ns t by the 
time (see eq. (3.35)) 

_n li i 11 _5_3£ 1 

Qinst ~ &s,inst(ichange) => Change ~ OL q^n h ~ £2 a 2 2 Q . (3.40) 

This occurs if 

l+c l+3c 

^iso VPmax j ^ ^-change ^ ^largc angle ^ Q 5 < 6 < Of 7 . ^o.4iJ 

This is the Region 3 in Fig. 10. 

In this regime the evolution of p max may cease to be controlled by saturation before 

l + c 

^change- In particular if e > a is , elastic scattering among the soft particles will 
start to push p max to higher scales (reducing T* as in the isotropic case) before the 
new instabilities take over momentum diffusion. However, q remains dominated by 
particles at the scale ki so , so this can have no effect on £ cn ange- 

• The third possibility is that soft particles never grow to dominate gj nst and the large 
angle changes erase the anisotropy, which occurs if 

l+3c 

e > a,— . (3.42) 
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Figure 10. For weak anisotropy, the d, c plane (representing the level of anisotropy and of over- or 
underoccupancy) divides into four regions. In region 1, plasma instabilities induce particle joining 
processes which erase the anisotropy. In regions 2 and 3, particle splitting creates soft particles 
which are anisotropic and generate new plasma instabilities. In region 2 this occurs because of 
particles which have saturated occupancies. In region 3 it occurs later, when the saturated particles 
have become isotropic and higher energy, still-anisotropic particles arc responsible for the new 
anisotropy. In region 4, the particles isotropize due to instability-induced angle change. Below 
the blue triangle, angle change is dominated by ordinary clastic scattering rather than by plasma 
instabilities. 



This is the Region 4 in Fig. 10. In this case the anisotropy shrinks appreciably at 

l+3c 

^change ~ *iarge angle given in eq. (3.19). When e falls below a~^~ the system will enter 
region 3. 

These conclusions are displayed in Fig. 10. 

All we have computed above is the time scale before some new physics significantly 
changes the dynamics; we have not computed the final equilibration time t cq . We will 
postpone doing so until the end of the next section. 

As a final aside, note that the total screening of soft particles actually dominates the 
screening arising from hard particles before £i ar ge angle if there are enough particles at the 
scale p max at that time. The total screening from such particles is m 2 ~ ap max T» ~ 
aPmaxQ- This is to be compared with m 2 ~ a 1 ~ c Q 2 from hard particles. Using eq. (3.28) 

l+9c 

and eq. (3.19), this occurs if e < a s . However this effect is not important because even 
if soft particles dominate screening, if they are isotropic they have no influence on which 
modes are unstable to plasma instabilities and how large the plasma instabilities grow. 
Therefore this large contribution to m 2 does not actually influence the physics of plasma 
instabilities and can be ignored. 



4 Large anisotropy 

A nonabelian plasma with a violent creation mechanism, or under strong anisotropic ex- 
pansion or compression, can be very far from equilibrium in a strongly anisotropic fashion. 
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In particular, we expect very strong anisotropy at early times in heavy ion collisions, at 
least in the theoretically clean case of extremely high energy nuclei (so that the energy den- 
sity after the collision represents a scale large compared to Aqcd and asymptotic freedom 
ensures perturbative couplings) [12]. Therefore we will consider very anisotropic systems, 
which may also possess typical occupancies very far from 1. 

Since we cannot study everything, we will concentrate on two of the simplest possibil- 
ities; an initial momentum distribution with a single characteristic scale Q and a strongly 
oblate (planar) momentum distribution, and a similar distribution but with a very prolate 
(linear) momentum distribution. We will try to determine the evolution of such a system, 
starting with the growth of the plasma instabilities, and continuing until some other piece 
of physics takes over as the dominant mechanism driving thermalization. We will not follow 
the system's evolution all the way to thermalization for every case we introduce, but at 
the end we will give an estimate of the final equilibration time scale t eq . 

4.1 Instabilities of an oblate distribution 

Consider first a distribution of excitations with a characteristic momentum p ~ Q, but 
with a very oblate angular distribution. Namely, the occupancy is f(p) ~ a~ c if p ^ Q and 
Pz ^ $Q, but f(p) is small outside this range. That is, the excitations are concentrated in an 
angular range of width ~ 5 about thep z = plane (so in spherical coordinates, | cos(0)| <, 5 
for most of the excitations). We take the angular range 5 to be small, characterizing it 
when convenient as <5 ~ a . Thus, d < (last section) represents distributions which are 
nearly isotropic, while d > represents distributions which are very far from isotropic. 
The energy density of this system is 

e ~ f d 3 p pf(p) ~ Q 4 Sa~ c (4.1) 

where the factor 5 is because the phase space volume which is occupied is only Q s 5, not Q s . 
If the energy density is larger than a~ 1 Q 4 then there will be Nielsen-Olesen instabilities, 
see subsection 2.3. This occurs if 5a~ c > a -1 or c — d > 1. We have already discussed the 
Nielsen-Olesen instability so we will not consider this region further. 
The screening scale is 

2 „ f dp ,, . ,J.-cxr& „l-c+d/ 



m 2 ~ a / — /(p) ~ a l ~ c 5Q 2 ~ a l ~ c+d Q 2 (4.2) 

and magnetic modes are unstable even in a larger range of k z . According to Ref. [25], the 
range of wave numbers which exhibit plasma instabilities is k such that k x ,k y < m and 
k z < 5'^ 1 m, and the growth rate is 7 ~ m. We can quickly re-derive these results by turning 
again to eq. (3.6) and eq. (3.8). If k x /k z < 5 and k y /k z < 5 then v ■ k < 5k for most of the 
excitations in the plasma. If 7 > v ■ k over most of this range then the contribution from 
the second term in eq. (3.8) is negative and dominant, and IIjj ~ m 2 k 2 /^/ 2 . The largest 
this can be, given 7 > v ■ k ~ 6k, is IL^ ~ m 2 5~ 2 . Substituting into eq. (3.6), the largest 
value of k which can be unstable is k 2 ~ II ~ 5~ 2 m 2 or k ~ 5~ l m. Since we required 
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Figure 11. Illustration of how a magnetic field can vary rapidly in the z direction, and yet particles 
will have long free paths in same-sign regions of magnetic field because they propagate in a narrow 
range of angles. 

k x ,k y & 5k and 7 ~ 5k in finding this largest unstable k value, we find that the transverse 
components of k must be ^ m, and the growth rate is 7 ^ m. 

Let us develop a clearer understanding of this result. We saw previously that the 
physical meaning of the scale m is that excitations must remain coherently within a region 
of magnetic field of stable sign and magnitude for a time scale 1/m if the current is to 
grow big enough to modify field evolution. It is immediately clear that 7 < m, since 7 
is directly an inverse time scale for the magnetic field to change. Now k determines the 
spatial nonuniformity of the magnetic field, and v ■ k tells how fast a particle explores this 
nonuniformity. If the field varies only in the z direction but a particle moves purely in 
the xy plane, it will remain in a region of coherent magnetic field forever. Since particles 
propagate almost purely in the xy plane, the k x and k y components must be ^ m so that 
in-plane propagation remains in a region of coherent magnetic field. Since the z component 
of an excitation's propagation is small, the field can vary faster in the z direction without 
the particle feeling different signs of B field. This is illustrated in figure 11. 

There is an assumption we have built into the above discussion. We have assumed 
that the excitations which give rise to these plasma instabilities can be treated as point- 
like entities. This only makes sense if they can be described using wave packets with 
spatial extent Az < l/k z , requiring momenta Ap z > k z . Clearly Ap z cannot be larger 
than 5Q, the width of the p z distribution of our excitations. Therefore our treatment is 
only self-consistent if 5~ 1 m < 5Q. Using eq. (4.2), this requires 
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m < 5 Q and m ~ 5*a 2 Q => 5 d > a or c < 1 — 3d. (4.3) 

Above this line, we must start from eq. (3.7) and re-derive what modes are unstable and 
how fast they grow. We will return to the case c > 1 — 3d at the end of this section. In the 
remainder of this section we will assume c < 1 — 3d, so the hard- loop treatment is valid. 

As in the previous section, plasma instabilities give rise to large infrared magnetic 
fields, which can dominate the dynamics by deflecting the hard excitations, with a char- 
acteristic momentum diffusion strength Qi ns t- We want to determine gi ns t, which requires 
understanding how large the instabilities grow. The plasma instabilities will grow until 
either 
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1. they cause some physics which changes the distribution of excitations or otherwise 
changes the dynamics in an essential way, or 

2. they get large enough that nonlinear effects limit their growth. 

Case 2 occurs first over almost all the parameter space we explore, so we will consider it in 
more detail. There will be a narrow region where case 1 actually occurs first, and we will 
handle it when it arises. 

The growth of instabilities will be cut off when any mechanism either interferes with 
the process which causes growth, or provides an efficient mechanism to remove energy 
from the unstable modes. Color randomization is a mechanism to cut off the growth 
mechanism, just as in the weak-anisotropy case. Refer again to Fig. 9. When the Wilson 
line of a propagating particle randomizes its phase over a length scale 1/v ■ k ~ 1/m, then 
color randomization occurs and instability growth is shut off. This occurs if the typical 
amplitude for a gauge field (associated with the plasma-unstable B fields) is A ~ m/g. 
The magnetic field is B ~ V x A ~ h mst A ~ m 2 5~ 1 g~ 1 , corresponding to 

B 2 ~a- l 5~ 2 m 4 . (4.4) 

Equivalently, we need a Wilson loop of size 1/m x 1/fc to contain an order-1 phase, requiring 
B ~ km/ g, returning the same estimate. Other estimates for the limiting size of the 
magnetic fields, such as the scale where nonlinear interactions between soft B fields become 
important, give the same estimate. This estimate is also supported by lattice studies [25]. 
The limiting size of the instabilities again corresponds to occupancies f(k) ~ a -1 for 
the unstable modes. To see this, note that the energy density in soft magnetic fields is 

/fcinat F~m 

dk z / d 2 k ± kf(k) ~ k 2 nst m 2 f{k) ~ m 4 5- 2 f{k) (4.5) 

requiring the typical unstable field occupancy to be f(k) ~ a~ l . 

These magnetic fields give rise to momentum diffusion, but now the momentum dif- 
fusion coefficient q is strongly direction dependent. For a particle moving at an angle 
9 > 5 with respect to the xy plane (that is, with p z /p = sin#), the magnetic field changes 
orientation on a length scale k\ ns \,/6, and so 

For generic angles this is m 3 /5, while for very planar directions it is m 3 /5 2 . The latter case 
is relevant for the typical excitations with p ~ Q, p z ~ 5Q. So for the typical excitations 
we have 

<7inst(<5) - <T 2 m 3 ~ rsai^Q 3 . (4.7) 

The situation becomes still more complicated if we consider relatively low-momentum 
excitations, with p < m/5 2 . Specifically, for excitations with p z < m/5, the extent of the 
wave packet is larger than the inhomogeneity scale of the magnetic field and the particle- 
approximation (WKB) treatment of propagation, implicit above, breaks down. Therefore 
we must consider p9 < m/5 as a special case. 
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The key to understanding the response of such a low-momentum particle is to re- 
member that the momentum transfer has to be a multiple of the wave-number of the 
nonabelian field. But if a particle of momentum p <C m/5 2 picks up a single kick of mo- 
mentum k z ~ m/5, then its angle will instantly change to be 9 ~ k z /p ~ m/{p5) S> 5. In 
this case, one should estimate the rate of momentum accumulation using this larger value 
of 9. Using this estimate, we find q ~ m 3 /59 but 9 ~ m/pS, so q ~ pm 2 . This estimate is 
valid for p > m/5. For p < m/5, there are /c-values which the momentum-p particle simply 
cannot pick up, since the magnetic field is nearly static and does not change the particle 
energy. Hence the largest k which can be picked up is k ~ p. Again, such a k would cause 
large-angle change to the particle, so the coherence time is only 1/k. This again gives the 
estimate 
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q ~ pm , when pm < -— . (4.8) 

59 

4.2 Momentum broadening 

Now we are ready to determine the range of momenta where plasma instabilities are more 
important than ordinary elastic scattering. The rate of ordinary elastic scattering is 



talc ~ a 2 [ d*kf[l+f] ~ a 2 - c 5Q 3 [l+f] ~ ( ^'Z^ ° H ' (4.9) 

J a 5Q c < . 



Later we will often find it convenient to work only in terms of the variables a,5,m,Q. Using 
that m 2 ~ a J f(p) cfip ~ af5Q s , we can eliminate the factors of / in favor of factors of 
m 2 and write 

Elastic (c < 0) ~am 2 Q. (4.10) 

Now we need to compare Elastic with <ft nst to determine where plasma instabilities dominate 
the dynamics. Within the range we are considering, d > and c < 1 — 3d, we find 
<?inst > ^elastic whenever c > 0, while for c < the instability dominates provided that 



ftnst > ^elastic 

al-f -IQ 3 > a 2 ~ c+d Q 3 

c > -l-3d (4.11) 

or equivalently d > - „~ c . Outside this region, elastic scattering is more efficient than 
plasma instabilities. Note that the region where elastic scattering wins out corresponds to 
systems with extremely low occupancy, c < — 1 or f(p) < a. We will not try to determine 
the physics in this region. 11 Also note that, for excitations at generic angles 9 ~ 1, the 
plasma instabilities give a smaller q, q ~ vr? /5 from eq. (4.6). For such particles, plasma 
instabilities dominate in the narrower region c > —1 — d. But this fact will turn out not 
to play a role in our discussion below. 



11 One might expect that the physics here will be analogous to isotropic, low occupancy systems studied 
in section 2.5. This may be true but it is not obvious, since low-momentum radiated daughters may be 
anisotropic and may generate important plasma instabilities. We leave this problem as an exercise to the 

reader. 
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Having identified the region where plasma instabilities dominate q and having found 
the relevant q, it remains to determine what physics first leads to a significant change to 
the dynamics. We will take as our criterion that some new physics causes q for our typical 
excitations to change appreciably. The simplest possibility is that the q we just found 
causes the typical excitations' angular distribution to get broader, increasing the value of 
5 and reducing the value of c. Since the mean squared change in p z is 

Ap 2 z (t)~qt, (4.12) 

this occurs when Ap 2 (t) > p 2 ~ 5 2 Q 2 , which requires a time 

ibroaden ~ ^ ~ ^ ~ ^a^^Q^ . (4.13) 

Self-consistency requires t > m _1 . We see that this is satisfied but is marginal for the 
boundary case c = 1 — 3d, see eq. (4.3). 

4.3 Splitting in highly anisotropic plasmas 

Just as for isotropic and weakly anisotropic systems, we also expect joining (for c > 0) 
and splitting (for c < 0) may be important. As for weak anisotropy, these processes will 
generically be in the LPM regime, so an analysis in terms of q and formation times is 
relevant. And the split daughters can again dominate the dynamics, either through elastic 
scattering or through new plasma instabilities. 12 

Assuming qi nst is the dominant momentum exchange process, and estimating the rate 
of LPM modified splitting or joining in the same way we did in previous sections, the 
formation time to split or join to the momentum scale p ^ Q is 



o . p 

tformiP) ~ \ 

Qmst 



( PS 2 ^ X-2 

^r P> z m, 



(4.14) 
— r v < 5~ 2 m . 



The different value when p < 5~ 2 m happens both because the momentum broadening 
rate is smaller for such soft momenta, and because these emissions are not actually LPM 
suppressed; 1/m is the mean time between scattering events so the emission rate will be 
the same as if we assumed each scattering [every time 1/m] provides the opportunity for 
an emission. The rate for a hard excitation to emit an excitation of momentum p is 

~^f p - at^ m [l+f{Q)} = r split (p) . (4.15) 

Here [1+/(Q)] is a Bose stimulation factor for the case that the "hard" p ~ Q excita- 
tions are at high occupancy. In this case, the final-state hard particle is Bose enhanced. 



12 In [24, 27] the authors argued based on numerical simulations of plasma instabilities that k <; m exci- 
tations are also produced by re-scattering of the unstable fields (the "nonabelian cascade"). While true, we 
find that this production method is less efficient than bremsstrahlung emission, a process which the lattice 
techniques employed in [24] do not incorporate. Therefore we will neglect the nonabelian cascade in the 
remainder of this paper. 
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(Bose enhancement for the emitted, soft particle is canceled by absorption, as discussed in 
subsection 2.5, see eq. (2.43).) 

The time scale on which a typical hard excitation undergoes hard splitting or joining 
is 

tspiit(Q) ~ a-V> c i form (Q) 3 V ~ a"+(i.4h3Q-i . (4.16) 

Here (0, c) refers to the low-occupancy or high-occupancy case respectively. This time scale 
is shorter than ^broaden if 

5*a 4^U'4j c < 52Q- 2^ W 

(5-S <ai + (t'T) c . (4.17) 

If 1 > c > 0, in which case we use the — c/4 case, this expression is never satisfied and 
broadening always happens before hard joining. Splitting will also play no important role 
when c > 0; just as in subsection 3.4, soft occupancies will become an f(p) ~ a~ c Q/p tail, 
which however never dominates screening, scattering, or plasma instabilities. Therefore if 
c > then the first new physics is always momentum broadening at time ^broaden- 

However if c < then in the region 5 > a( _1_3c " 9 , or c < —3d— 1/3, hard splitting will 
occur before the plasma instabilities can broaden the hard particle distribution. Therefore 
splitting will play a role in a region at least this large. The actual region will be larger, 
because soft p <C Q split-off daughters can dominate the dynamics before hard splitting 
becomes important. 

4.4 Scales h mo and k 8p ut 

The analysis of what physics takes over from the "ordinary" instabilities will be a little 
intricate. In particular, one complication is that, as we have seen, the momentum broad- 
ening rate <ftnst(0) is angle-dependent. Since we will often be concerned with whether or 
not the radiated daughters are isotropic, this is a complicating effect. We summarize how 
the daughters fill in small momentum scales, including this complication but neglecting 
splitting of daughters, in Appendix A. 

However, in the end we are always interested in what is happening "at the last minute" 
when some other piece of physics comes to dominate q(9 ~ 5). In every case we will find 
that this new piece of physics has a q which is not parametrically anisotropic. That is, 
q will be anisotropic by an 0(1) amount, but not by a power of 5 or a. Therefore, at 
the time scale when a new q is coming to dominate, we may treat q as approximately 
isotropic. Similarly, in eq. (4.14) we can ignore the p < 5~ 2 m case, which occurred because 
of the anisotropy of q, and always take tf OTm to comply to the first expression. It is most 
important to get this last time scale correct. If we also treat q as isotropic at all times, 
while we will mis-represent physics at early times, we always get the final stage, the time 
scale and physical mechanism which replace the original plasma instabilities correct. 

To understand how radiated daughters can come to dominate q, we need to study the 
time development of the population of such daughters. We will distinguish two important 
scales. The first is the scale ki SO (t), introduced in eq. (3.29). This is the time-dependent 
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momentum scale such that particles at this scale have had just enough time for transverse 
momentum diffusion to drive them to an isotropic distribution at the order- 1 level. It is 
determined by 

kf so ~ <?inst* -> hso ~ <5 _1 m2i2 ~ a t t*Q* . (4.18) 

Excitations at the scale k lso can grow to dominate q by providing a new source for plasma 
instabilities. Unlike in the case for weak anisotropy, the scale p max lies below ki so (t) for all 
t > m . The scale p max can only play a role by dominating q through elastic scattering. 
However, if elastic scattering ever provides ^elastic ~ <Zinst, then all excitations with k <^ ki SO 
thermalize and ^elastic is controlled by the scale /cj so . Later we will check whether elastic 
scattering from the scale /c; so can dominate, and we find that this is not the case in the 
regime of interest. 

The other important scale is k sp n t (t), defined in eq. (2.71): r s piit ( ^spiit ) * ~ 1- As 
discussed in subsection 2.5.3, a daughter with k < /c S piit i n turn has time to re-split and 
fragment completely. Therefore, all particles with k <> k sp m redistribute into a thermal 
bath or thermal tail. Applying eq. (4.14) and remembering c < so [1+/(Q)] ~ 1, we find 



Split 

3 
771,2 



^^"piit-^Vm-a- 1 ^^ => k SV)Vlt ~5- 2 a 2 mh 2 -oT^Qh 2 . (4.19) 



Both k\ so and £; S piit increase with time; but fc S piit rises faster, as t 2 rather than as t*. 
The scales cross at a momentum fci SOS piit an d time £i SO spiit such that 

"■split V^isosplitj "'isov'-isosplitj = ""isosplit 

3 i 

<T a 2 m 3 t? osp iit ~ ^"^Lpiit 

^Lspiit ~ 5a ~ m ~~ 2 

iisosplit ~ 5la-lm- 1 ~ a -T+i c +\ d Q- 1 , (4.20) 

^isospiit ~ ^~ 2 a 2 ^ 3 4os P nt ~ S'sa'sm ~ a e Q . (4.21) 

At all times, modes with k > fcisosplit split before q drives them to isotropy, while modes 
with k < fcisospiit become isotropic before splitting. 

4.5 t < iisospiit : new instabilities 

Before the time iisosplit , &iso > ^spiit, and we can effectively ignore re-splittings of daughters. 
The daughter particles are generically anisotropic and can give rise to "new" plasma insta- 
bilities. We will now show that the most important scale for any "new" plasma instabilities 
is the scale fe; so . Due to the nonabelian LPM effect, the number of particles varies with 
scale as f{k)d?k ~ k^ 1 ' 2 dk/k which is IR dominated; the contribution to screening scales 
as f(k)d 3 k/k ~ k~ 3 ' 2 dk/k and is more IR dominated. Define the contribution to screening 
from momenta around k to be m 2 (k); we see m 2 (k) ex k~ 3 ' 2 . The level of anisotropy also 
varies with scale. For k > ki so the particles are strongly anisotropic but are few in number. 
Since the mean squared transverse momentum they can accumulate is kf , they will lie in 
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Figure 12. Scales k\ so and /c S piit intersect at time ii SOS piit at momentum scale fci SO spiit- All modes 
above fci SOS piit (e.g. the scale Ua in the plot) split before they isotropize, while modes with k < /cisospiit 
(such as the scale ks) become isotropic before splitting. 



an angular range 6(k) ~ ki so /k. These daughters will give rise to new instabilities with a 
q of order 

q(k > k iso ) ~ 0- 2 m 3 {k) oc {k/k iso ) 2 k- g / 4 oc k~ 1/4 (4.22) 

which is dominated by small k. Meanwhile, for k < &i so , the excitations are nearly isotropic. 
Only the fraction which were introduced too recently to isotropize are anisotropic, and 
this fraction is of order {k/k\ so ) 2 . Therefore, for k < k{ so , we have weak anisotropy with 
e ~ (k/ki so ) 2 . Then we find 



q(k < k i80 ) ~ e2m 3 (k) ~ (k/kisofk * oc fc< 



(4.23) 



which is large k dominated. Hence the most important scale is k[ so , where the anisotropy 
is order- 1. 

To determine the value of q(h so (t)) = q S) i ns t due to plasma instabilities at the scale 
hso(t), we need to determine m 3 (ki so ). We have that t 2 oim (hso) ~ hsoQ" 1 - But q ~ k 2 so t~ l 
so t 2 orm ~ t/ki so ~ 5t2m~2 . The number of daughters produced is n(ki so ) ~ ^hard^i/iform- 
And m 2 ~ an^^/Q so an^d ~ m 2 Q. Therefore 



n \k[sc 



^form("'iso, 



m 2 Qt 2 , i i 



(4.24) 



and m 2 (k iso ) ~ an(/c iso )//c iso , so 
m 2 (k iso ) ~ am 2 fe iso 2 Qt2 ~ a^m^Qi* =>. g sinst ~ m 3 (ki SO 

The strength of plasma instabilities from daughters grows with time as t 



3 3 15 3 3 

a2 i5*m s Q 2 £s 



(4.25) 
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These "new" plasma instabilities will catch up with q{ nst ~ 5 2 m 3 at the time scale 



3 3 15 „ 3 



a5(54m '*Q~ 2 t* mst ~ 5 2 m 3 
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2 



t s ,inst ~ a~ 4 <5~m 3 Q~ 4 

~ a -f-f c -f ^Q- 1 . (4.26) 

Of course, we assumed £ S) i ns t < fisospiiti so this answer is only self-consistent if 

5 3 „ 35 j „ i 11 i c i d „ i , 1 , _. 

Wnst < tisosplit =* a-^-s^-^Q^ 1 < a-ir+s+eQ- 1 => c <_3d--. (4.27) 

Below this line in the c, d plane, new instabilities from split daughters dominate q before 
p ~ Q particle angle-change becomes important. This happens to be the same as the line 
we found earlier, when we determined where hard splitting happens before ^broaden- 

We should also check that the time scale we found in eq. (4.26) is longer than the 
time scale 1/m. After all, the instabilities grow on this time scale, and it is the shortest 
formation time of any daughter induced by scattering in the instabilities. Therefore, if we 
find t Sj i nst < 1/m, all of our calculations are inconsistent. The condition is 

-1 _4 r z22 3^-4 -1 -2-2c-^-d -, r — 3(1 + c) 

(4.28) 
which is a line starting at d = 0, c = —1 and rising with slope —3/8. Between this line and 
the line d = — (1 + c)/3, daughters become important and control q via plasma instabilities 
in a time scale ~ 1/m, possibly before the plasma instabilities from the original hard 
particles have finished growing. We have not completely worked out the physics in this 
region, and we suspect that the physics in the region immediately below the d = — (l + c)/3 
line may also be nontrivial. However we will be lazy and leave the complete understanding 
of this region for future work. 

4.6 Elastic scattering versus new instabilities 

Let us return to the region with d < —(1 + 3c)/9 but d > —3(1 + c)/8. We just found that 
new plasma instabilities, associated with the scale ki BO , dominate q before ii SO spiit i n this 
region. However there is one more thing we should check. The modes at and below the 
scale /ci so could in principle also dominate q via elastic scattering. Before continuing, we 
should pause to make sure that this does not occur. 

Suppose to the contrary that ^clastic ^ q{hso) which we just computed. Then /ci so due 
to elastic scattering will be at least as large as that from plasma instabilities. But as 
discussed in subsection 2.5, elastic scattering will fully thermalize all k < ki so , either into a 
thermal bath or an f(p) oc T*/p type tail. We will consider separately the case f(h so ) S> 1, 
in which case we get a T*/p type tail, and the case f(ki SO ) <C 1, in which case we get a 
thermal bath. 
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Begin with the case f(ki so ) S> 1, and compare the sizes of g e iastic with q s ^ ms t: 

/fciso 
/[!+/] *k whereas * S)inst ~ ^ 

Elastic ~ « 2 / 2 fcL whereas q s ; mst ~ (afkf so )i ~ a^/ 1 fcf so . (4.29) 

We see that Elastic < <?s,inst at time i Siinst provided that /(fcjso(*a,inst)) < a" 1 - We can also 
quickly check that ^elastic (*) was not larger at some earlier time; we have fci so ~ t?, and a 

_7 3 

quick calculation shows f(ki so ) ~ A; iso 2 t ~ t~*, and therefore g e iastic ~ t°. 

On the other hand, if f(k iso , t s ,i nst ) < 1 and again assuming that <? e i astic > q s ,i m t, 
then elastic scattering re-organizes the momentum distribution of k < ki SO particles into 
a thermal distribution. The temperature is set by the energy density, T 4 ~ f(ki so )kf so . 
Meanwhile, particles with k just above ki SO have not "collapsed" and still provide q ~ q s ,mst 
via plasma instabilities. Making the comparison again, we find 

Elastic ~ a 2 T 3 whereas q s , inst ~ m 3 
Elastic ~ a 2 (fkf so )4 whereas q s , inst ~ (a fkf so )z 
^elastic ~ a 2 / 3 kf so whereas q S: - mst ~ 02/2 kf so . (4.30) 



In this case, q St mst is larger provided f(ki SO (t s ,imt)) > a 3 - In this case T grows with time, 
and so Elastic wa s smaller at all earlier times. Combining the two cases, ^ s ,inst dominates 
provided that a >> / 3> as . 

Now let us check the actual value of /(A;i so ,t s ,inst)- Using / ~ n(ki so ) / kf so and combin- 
ing eq. (4.18), eq. (4.24), and eq. (4.26), we find 

K iso 

(4.31) 
This is to be applied in the region — (1 + c)/3 < d < —(1 + 3c) /9. In this region, / varies 
from a -1 , at d = —(1 + c)/3, to as, at d = —(1 + 3c)/9. Therefore we self-consistently 
find that new instabilities dominate elastic scattering when they come to dominate overall 
q, throughout this region. 

4.7 t > tisospiit : new thermal bath 

Now consider the region (1 — c)/3 > d > — (1 + 3c)/9. In this case, at the time iisosplit, the 
total q is still dominated by gi nst arising from instabilities induced by the p ~ Q excitations. 
We want to determine whether split daughters can change this at times iisosplit < t < 
^broaden- This range of momenta only exists if £ is0 split < ^broaden, which, using eq. (4.13) 
and eq. (4.20), requires d < — (1 + 3c)/7. 

Below the scale fc sp iit, excitations undergo split-join processes on a time scale shorter 
than the age of the system. This allows the excitations with k < k 8p m to equilibrate into a 
nearly-thermal bath with temperature T, in analogy with what we found in Sub-subsection 
2.5.3. Let us find T(t) as a function of t. The thermal bath is made by breaking down all 
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daughters produced with k < kgpin. As in subsection 2.5, most of the energy comes from 
the scale fc sp ii t ; the number of daughters is ~ ^hard; so 



-U _ ~2r - x -- 1 - 1 



£ ~ fc split n hard ~ a k spht mQ, T ~ £4 ~ a 4 fc s 4 plit ra 2 Q 4 . (4.32) 

Substituting in eq. (4.19), 

T ~ aU-^rniQh? . (4.33) 

This grows more slowly with time than fc sp iit- We can also explicitly check that at the 
earliest time we are considering, namely ii SOS piit; 

Tfeospiit) ~ a~T5<HmtQ3 ~ a" 1_ 24 +M Q (4.34) 

whereas, from eq. (4.21), 

_2 _2 -4-12c-4d 

Msosplit ~ a 3^ 3m ~a 24 Q. (4.35) 



Since we are considering a region where 3c > — 1 — 9d, we find T(£ isosp i it ) <C Msospiit ■ Hence 
T <C A; sp ii t at all subsequent times. We also find that T <C fci so at all times, since fei so = fcgplit 
at time £i SO spiit> an d T and ki so scale with time in the same way, as £2 . 

Most of the particles in the system are in the thermal bath. To see this, note that each 
particle at energy k sp nt turns into ~ (k sp n t /T) S> 1 thermal bath particles. Since ~ n na rd 
particles of energy k sp n t get made, the bath particles outnumber hard particles, k ~ fc S piit 
particles, and (as we will soon see) particles with k sp nt S> k S> T. Hence we expect that 
elastic scattering will be dominated by the thermal bath. The value of q due to elastic 
scattering with thermal bath particles is 

<7 s ,ciastic ~ a 2 T 3 ~ a^d-^m^Qiti . (4.36) 

In addition, the thermal bath is generally not fully isotropic, and neither are any of 
the excitations at scales > T. Therefore there can also be plasma instabilities due to the 
anisotropy of excitations at any scale between T and fc S piit- First we need to estimate 
occupancies in this range. We know n(k sp y lt ) ~ "-hard- What about scales below fc sp iit? 
Below the scale fc S piit j the flux of energy density through a logarithmic "bin" of momentum k 
is ~ & S plit n hardA- The energy density equals this energy "flux" times the mean lifetime of a 
particle to undergo democratic splitting, 13 which is £ S plit(fc) °c & 5 ( see ec L- (4-14), eq. (4.15)). 
Since the energy density is kn(k), we find n(k) oc fc~2, specifically n(k) ~ rihard(&split/&) 5 - 
The strength of screening from these particles is m 2 (k) ~ a~ l n(k)/k oc /c" . 

Recall that, for k > fci SO spiit, particles split before their directions are randomized. Since 
fcsplit > fcisospiitj the largest k values in the cascade will be highly anisotropic. Specifically, 
in the time t sp \it(k) during which a particle resides at momentum k before splitting further, 
it picks up a mean-squared momentum qt S plit(k). This induces an angular range of 

d 2 {k) „ ^! „ «te K k -l . (4.37) 



13 By democratic splitting we mean a splitting process where the daughters have comparable energy. This 
is the main process which removes particles from some fc-bin. 



43 



There is also a contribution to 6 2 from the angular distribution of the parent, but since 
we find 2 increases with smaller k, the parent-distribution is subdominant and we can use 
the above estimate. The strength of plasma instabilities from a scale k 3> fci sosp ii t is 

?s,mst(fc) ~ 9- 2 (k)m 3 (k) oc k%k~* ~ k'i (4.38) 

which is small k dominated. 

If T < feisosplit) then there are also nearly-isotropic scales between fc sp iit and T. The time 
scale for a particle of momentum k < fei SO spiit to be direction-randomized is t{ so (k) ~ k 2 /q. 
A fraction t- mo (k) /£ S plit(&) of particles have not had time to randomize in direction, giving 
rise to a residual (weak) anisotropy of size 

<k) ~ — ^— oc ki . (4.39) 

g*s P iit(«) 

The plasma instabilities arising from such excitations are of order 

q s ,wst(k < k isospht ) ~ (em 2 ) 2 oc A; , (4.40) 

essentially flat in k. Therefore it is fair to say that, up to logs, the size of new plasma 
instabilities is given by the contribution from the softest scale in the game, which is T 
itself. 14 

At the scale T, the time it takes for an excitation to randomize its direction is 
t; so (r) ~ T 2 /q. Any excitation which arrived longer than ti so (T) ago has had its di- 
rection randomized; recent additions to the bath have not and are anisotropic. Therefore 
the degree of anisotropy of the bath is 

tiso(T) T 2 a^8~ l m^Q^t i .-, _i^i i+c+zd 

em ~ -^- ~ ^— ~ -r-^ — 5 ~ a 25m 2Q2 ~« 4 (4.41) 

w t qt 5" 2 mH ^ y ' 

which turns out to be t independent. 

If we consider the thermal bath as a weakly-anisotropic system, in the spirit of section 3, 
then it has c = 0. In this case, we saw there that instabilities dominate q from the thermal 
bath if e > as. Otherwise elastic scatterings dominate. Therefore plasma instabilities are 
more important provided 

1 i+c+3d 1 1 — 9d . , ,„. 

e > «3 => a 4 > a s => c < — - — . (4.42) 

In this region, the value of q from T-bath plasma instabilities is 

9s,inst(*) ~ e 5 "^s ~ eaa^T 3 ~ a 3 m 3 Q2t2 . (4.43) 



14 If T < fci SOB piit, the argument is less clear. In this case, the excitations may isotropize as they cascade 
from the scale fci sosp ij t to the scale T, "landing" on the thermal bath already isotropic. However, in this 
case the q of plasma instabilities from the 0(1) anisotropy at the scale fci SO spiit turns out to be the same as 
the q we find for the thermal bath below, so our final answer remains correct. 
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Now we just need to check whether q from eq. (4.36) or eq. (4.43) can exceed gi nst ~ 
5~ 2 m 3 before the time t = ibroadcn- New instabilities catch up with qi nst at time £ Si i ns t, 
defined as 

Q's,instl,' : s,i ns tJ ~ 9inst 

a 3 m 3 Qhl inst ~ 5~ 2 m 3 

t B ,u*t ~ a~ 2 6-iQ- 1 . (4.44) 

Comparing with ^broaden from eq. (4.13), we find 

£s,i ns t < ^. ^broaden II 

a^d'^Q- 1 < 5 4 m- 3 Q 2 

o Ad _ i — 3+3c+5d -i 

> 3 + 9c + 23d (4.45) 

a line starting at c = —1/3 and with slope —9/23. Therefore, new plasma instabilities from 
the scale T* come to dominate the dynamics in the region with 9d > — 1 — 3c, 23d < —3 — 9c, 
and 9d < 1 — 3c. 

For 9d > 1 — 3c, the thermal bath is too isotropic to have important plasma instabilities 
but (/^elastic can grow to dominate. This happens at time 



^s, clastic v^s,clast j ~ <Zinst 



11 _ 3 15 3 J 



aid 5m4Q4tJ ^ 6 2 m 3 



s.clast 

11 _— 1 -1 _ -i 



^,clast ~ a e S 3 m 2 Q 2 . (4.46) 

Again comparing with tbroaden, we find 

£s,clast ^ ''broaden U 

-11 -1 -1 -1 A _o 9 

a e 5 3 m 2 Q 2 < (5 m d Q z 

-25+3c-7d _-. -18-j-18c+30d _-, 

a 12 Q < a 12 Q 

> 7+15c + 37d. (4.47) 

Therefore elastic scattering from a thermal bath dominates in a region satisfying 37d < 
—7 — 15c but 9d > 1 — 3c. This is a triangle in the c, d plane with its apex at c = —25/6, d = 
3/2. 

This completes our study of the region where plasma instabilities are important and 
are well described by hard loops. The regions we have found and their properties are 
summarized in figure 13. The last region studied, where elastic scattering from a thermal 
bath comes to dominate, lies off the edge of the figure. 

4.8 Highly anisotropic oblate distribution 

Now we return to the case c > 1 — 3d. Recall that, in this region, the hard- loop approx- 
imation returned a range of unstable modes which is broader than 5Q. This indicated 
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an inconsistency in the hard-loop approximation. We also found for the marginal case 
c = 1 — 3d that the time scale ^broaden f° r the momentum distribution to broaden was the 
same as the 1/m time scale for the instabilities to grow. Let us try to understand the 
physics in this regime. 

We claim that the dominant physics remains a Weibel-type instability, with large 
magnetic fields deflecting the "hard" p ~ Q excitations and broadening the p± distribution 
of their momenta. However, the range of plasma-unstable modes will be different than 
what we found using the hard-loop approximation; and they will grow large enough to 
broaden the p ~ Q excitations' angular distribution before they grow large enough for 
nonlinear (nonabelian) interactions to cut off the instability growth. That is, the growth of 
instabilities will be regulated by the angular broadening of the hard excitation distribution, 
rather than nonabelian interactions between unstable modes. 

We will guess that the unstable modes still have k z S> k± ~ m but k z <C Q- But we 
now assume k z > p z ~ 9Q in contrast to the previous section. The one-loop self-energy is 
still given by eq. (3.7), but now p • k ~ Qk° + Qkj_ + p z k z need not be large compared to 
k 2 ~ k 2 . In order to find the range of unstable k z available, take k° and k± to zero; then 
p ■ k ~ Pzk z - But we expect k z > p z , so k 2 ~ k 2 is larger than p ■ k. So to find the range 
of unstable k z , drop the p ■ k terms in favor of k 2 terms. In this approximation, eq. (3.7) 
simplifies to 

n^(fc) ~ g 2 I d 4 p^S(p 2 )f(p) ~ %g 2 f d 4 P 5(p 2 )f(p) ~ %m 2 . (4.48) 

J K z K z J K z 

The propagator equation has instabilities when k 2 ~ n(fc), which we see requires 

Q 2 



k2 z~T2 m2 => k z ~^/mQ. (4.49) 

k z 

This is the range of unstable momenta. In the parametric regime considered, it is wider 
than p z ~ 5Q (so the treatment is self-consistent) but narrower than (5^ 1 r?i. 

What is the range of k± and the maximum growth rate 7? They are determined as 
the values which significantly reduce the value of W J,U estimated above. This occurs when 
p ■ k ~ kg. Since p° ~ Q ~ pj_, this occurs whenever 7, k± ~ m. So the transverse 
momentum range and growth rate of the unstable modes remain ~ m. 

To what amplitude do the plasma instabilities grow? As before, they grow until either 
nonlinear physics interferes with the growth mechanism, or the unstable fields modify the 
dynamics in some other way such that the calculation of the instability growth rate gets 
amended. In this region we will find that the latter occurs. Our previous estimate for 
the magnetic field strength where nonlinearities such as color randomization shut off the 
instability can be extended in a straightforward manner. The fields would still saturate 
with occupancy f(k) ~ a -1 , with a magnetic field strength of B 2 ~ a~ l k 2 m 2 ~ a~ 1 m?Q. 
The maximum value of q would then be m 2 Q. 

The unstable fields grow exponentially, so the momentum randomization will explore 
all values smaller than the above q before achieving this value of q. However, as soon as 
Qmst ^ S 2 mQ 2 , then the hard excitations undergo significant broadening, Ap± > p± ~ 5Q, 
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Figure 13. Cartoon of d,c plane indicating the regimes found for oblate distributions. Plasma 
instabilities are important in regions 1 to 7. In region 1, instabilities induce splitting which reduces 
anisotropy. In 2, instabilities induce splitting, generating anisotropic, occupancy-saturated parti- 
cles which cause stronger plasma instabilities. Region 3 is similar but the new particles are not 
occupancy-saturated. In regions 4 and 7, plasma instabilities directly lead to the broadening of the 
primary-particle distribution. In 4 they do this after amplitude-saturation, in region 7 they do so 
before their amplitude saturates and they are not well described by hard-loops. In region 5, splitting 
creates new anisotropic particles on a time scale of order the instability growth rate. In region 6, 
splitting generates a thermal bath, but it is incompletely isotropic and its plasma instabilities come 
to dominate dynamics. Region 8 is Nielsen-Olesen unstable. In region 9, elastic scatterings play the 
dominant role. Not shown is region 10, between regions 4,6 at very negative c and large d, where 
clastic scattering from a new thermal bath comes to dominate dynamics. The lower boundaries of 
regions 2,5 have not been completely studied and may require revision. 

in a 1/m time scale. In the region we are considering, m > 8 2 Q and so 5 2 mQ 2 < m 2 Q. 
Therefore the unstable modes reach a strength where they broaden the hard excitation 
distribution before they finish growing to become nonlinear. 

This broadening of the hard excitation distribution changes which modes are unstable, 
how large they can grow, and so forth. In fact, it moves where the system lies in the c, d 
plane. Therefore it constitutes an important change to the properties of the system. 

Hence we conclude that the physics of the region c > 1 — 3d but c < 1 + d is, that 

1 — c+d 

plasma instabilities grow in a time scale m ~ a 2 to become large enough that they 
broaden the distribution of hard excitations. The values of c, d then move along a line of 
fixed energy density, which is fixed c — d, towards smaller c, d, reaching c = 1 — 3d in a 
time scale t ~ 1/m up to logs. Thereafter they enter the c < 1 — 3d regime we have treated 
above. 

4.9 Final equilibration times 

In the previous subsections we have shown that plasma instabilities generally dominate the 
dynamics at early times, and we have determined what new phenomenon takes over from 
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them and after what time scale. The 10 regions we found can be broken into four broad 
categories: 

1. In region 9, elastic scattering dominates and the discussion of section 2 should be 
used. 

2. In region 8, N-0 instabilities almost instantly change the system to be one lying just 
off the boundary of region 8. The equilibration then proceeds as it would in that 
other region. 

3. In regions 1, 4, and 7 plasma instabilities cause the momentum distribution to become 
more isotropic. For d > this means c, d each decline with time; for d < only d 
declines with time. Therefore after some amount of time the system enters region 3, 
6, 9, or 10. 

4. In regions 2, 3, 5, 6, and 10, a bath of soft excitations develops and comes to dominate 
the dynamics. 

Since we already solved the first case, and since the second and third cases turn into either 
the first or fourth after some amount of time, we only need to address what happens in 
the last case, when a bath of soft excitations emerges to dominate the physics. This bath 
may start out nearly thermal (regions 6, 10) or far from thermal (regions 2, 3, 5). But in 
every case we have that the energy density e ~ a~ c Q 4 (if d < 0) or ~ a~ c+ Q 4 (if d > 0) 
is small compared to Q 4 , so the final temperature Tfi na i ~ e* <Q. 

In every case we expect the late evolution to be controlled by the physics of a bath of 
soft p ^ 7fi na i excitations, which steal energy from the hard p ~ Q modes and eventually 
become a thermal bath. Since the thermalization time of soft excitations is generally shorter 
than for hard excitations, one expects that the soft bath becomes nearly thermal before it 
finishes breaking down the p ~ Q excitations. It is safe to assume that, as in the isotropic 
case, the most time consuming stage is the final breakdown of the hard excitations when 
the bath temperature approaches T ~ Tfi na i. 

There are two possibilities. Near the end, the thermal bath can either be somewhat 
anisotropic, with e > Q3; or it can be more isotropic, e < as. In the former case the 
dominant physics will be plasma instabilities; in the latter case it will be elastic scattering. 
In either case, what we need to compute is the time scale at which the hard p ~ Q 
excitations fragment completely, that is, fc sp iit(£ C q) ~ Q- Using our previous results for 

"'split ) 

t eq ~ a^q'^Q^ . (4.50) 

It remains to self-consistently determine e and q. 

1 3 3 n 

First assume e > as so plasma instabilities matter. Then q ~ e2 « 2 ^fi na i- We also 
have that fci SO spiit ~ a~^~e 2 Tfi na i which is 3> Tfj na i; so the fragmenting daughters are nearly 
isotropic just above the scale Tfj na i. However, because q is caused by plasma instabilities, 
it is order- 1 anisotropic, and so the particles fragmenting into the thermal bath are order- 1 
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anisotropic. 15 The level of anisotropy of the thermal bath e is then set by the ratio of the 
time scale for direction randomization, ti so ~ T^ nal /q, and t eq - since a fraction ~ iisoAeq 
of the thermal bath is made up of excitations which landed less than ii so ago and have not 
isotropized. That is, 



tiso ^final? i 



3 _i „ 1 



"-cq a l q 2Q2 



and hence 

2 -13 „ 5 , 



-12 



e~c^(7fi nal /Q)r, t eq ~ a" Q^T fin 7 al . (4.52) 

1 2 

This result only makes sense if e > as, which we see is Tg na i > azQ. This is the region 

> c > —8/3 for d < and > c — d > —8/3 for d > 0. It also only makes sense if the 

resulting i cq > (x~ 2 T^ ngl . Otherwise, the p ~ Q excitations break down before the bath 

can completely thermalize with itself, and this final thermalization is the slowest process. 
2 
For Tg na i < a-iQ, thermalization takes so long that the thermal bath is very nearly 

isotropic and elastic scattering dominates. In this case q ~ a2 ^i na \j an d we find 

teq ~ a~ 2 Q^T^ . (4.53) 

This case is the same as the isotropic case from subsection 2.5. 

There are two weaknesses in this analysis. The first is that we assumed that there 

is no long time-scale "hang-up" on the way to the final thermalization we discuss. This 

seems safe; we checked explicitly that no such "hang-up" happens in isotropic systems, 

and anisotropic systems should generally have faster dynamics due to plasma instabilities. 

Also, all the £ c hangc time scales in Regions 1-7, 10 in table 1 are parametrically shorter 

2 
than t cq . The second weakness is that we assumed that, for Tg na i > asQ, the thermal bath 

will be anisotropic. Once it is anisotropic, the anisotropic arrival of daughters maintains 

this anisotropy, but it has to get that way to start with. However, in regions 2, 3, 5, and 6 

the soft bath does start out anisotropic; and region 10 (where it does not) is purely in the 

2 
region where Tfi na i < a-iQ. So this seems consistent. 

It is possible that, even if a system starts out with Elastic somewhat larger than q\ nS f, the 

thermal bath will still manage to become anisotropic. Therefore the boundaries between 

regions 2 and 9, and between regions 5 and 9, are somewhat in doubt. We will not try 

to resolve this issue here but leave it for future work, if physical systems which lie in this 

regime are found. 

4.10 Prolate distribution 

Now we will repeat the exercise of the last subsections, but for the case of a highly prolate 
anisotropic plasma, that is, one where f(p) ~ a~ c provided p z < Q and p x ,P y ^ SQ with 
5 ~ a d <C 1. 



15 This sounds like a contradiction; the particles are nearly isotropic as they fragment, but they land 
on the thermal bath in an anisotropic fashion. The reason this can be true is that, at every stage, the 
fragmentation happens in an anisotropic fashion dominated by particles near the p x ,p y plane, but they 
then isotropize before the next step in the fragmentation. 
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The number density of typical excitations is n ~ /<5 2 Q 3 , now with two powers of 5 
because momentum on two axes is more constrained. The energy density is 

e ~ nQ ~ a~ c 6 2 Q i (4.54) 

and we will assume cT c b 2 <C a" 1 , or c < 1 + 2d as otherwise there are Nielsen- Olesen 
instabilities. 

The screening scale is 

m 2 ~ anwiQ- 1 ~ <5V~ C Q 2 ~ a l - c+2d Q 2 . (4.55) 

Since the "hard" p ~ Q modes all travel nearly in the z direction, they remain coherently 
in a region of constant magnetic field provided that k z < m for the magnetic field; but 
both k x and k y may be larger by a factor of 5" 1 and the particles will remain in a region 
of coherent field. Hence we find that the range of unstable momenta and the growth rate 
are 

k z ~ m , k± ~ m5~ , 7 ~ m . (4.56) 

All unstable modes have polarization vectors in the z direction. The hard-loop approx- 
imation breaks down if k± for unstable modes exceeds the typical transverse momen- 
tum of a "hard" excitation, k± > 5Q, which is if m > 5 2 Q or c > 1 — 2d. The region 
1 + 2d > c > 1 — 2d is analogous to the case of subsection 4.8; in this region the growth 
of plasma instabilities is cut off before amplitude saturation because the hard momentum 
distribution gets broader. 

The unstable mode amplitudes saturate when, in traversing a distance 1/m along the 
z axis, a particle's color is rotated by an 0(1) angle, as previously discussed. This requires 
A z ~ m/g, B ~ kA ~ m 2 5~ l a~2. This is equivalent to the condition that a Wilson loop 
of extent A;" 1 by k~ l returns an order-1 phase. The magnetic energy density and q are 

B 2 ~ oT x n£b- 2 , q- mst ~ aB 2 /k z ~ m 3 <5" 2 . (4.57) 

Unlike the oblate or weak-anisotropy cases, saturation now occurs before the typical occu- 
pancy of an unstable mode is ~ a" 1 ; B 2 ~ fk^k z so / ~ a~ 1 5 for the unstable modes. 

The functional form of <7i nst is the same, when expressed in terms of m, 5, a, Q, as for 
the oblate case, eq. (4.7). The same will be true, for c < 0, of (feiastio as we now show. The 
momentum diffusion due to ordinary scatterings is 

^tic -o?\ d 3 P f[l+f] ~ ( « 2 ~W c > ' (4.58) 

J I am Q c < . 

For c > instabilities always dominate q and we can ignore g e iastic- For c < 0, eq. (4.58) 
has the same form as eq. (4.10). The region where plasma instabilities dominate is 

9inst > ^clastic 

a^-'W > a 2 ~ c 5 2 Q 3 

5- 1 > a 1 ^ , (4.59) 
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Figure 14. The same as figure 13 but for the case of a prolate anisotropic system. The meaning 
of each region is the same as in the previous figure; only the slopes of the lines separating regions 
have changed. 

or 2d > -1 - c. 

We now go quickly through the prolate versions of the arguments we gave in the 
previous sections for oblate systems. We will find that, when expressed in terms of m 2 , 5, a, 
and Q, the results from the oblate case carry over to the prolate case. Expressing in terms 
of c, d, the difference will just be the power of d appearing in to 2 , eq. (4.55). Therefore the 
distinct regions in the c, d plane are the same as before but the boundaries are found by 
replacing factors of d— c in the oblate case with 2d—c in the prolate case (since m 2 oc a 1+d ~ c 



in eq. (4.2) but m 2 



oc a 



l+2d-c : 



in eq. (4.55)). 



In particular, we find that eq. (4.13) for ^broaden, eq. (4.14) for tf orm and eq. (4.15) for 
tsplit, eq. (4.18) for k iso , eq. (4.19) for k spht , and eq. (4.20) and eq. (4.21) for * is0 split and 
Msospiit> all still hold if one reads the versions expressed in terms of m 2 , a, 5, Q. So do our 
arguments regarding the most important scales for new daughters. The temperature of 
the thermal bath of daughters, T*, is also still given by eq. (4.33), and the criteria for its 
dominance are determined by the same expressions when we use the to 2 , a, 5, Q variables. 

Therefore we find the same regions as for an oblate distribution. In terms of c, d the 
regions' boundaries are shifted, but in terms of d and the energy density (which depends on 
c — d for oblate and c — 2d for prolate distributions) they are the same. This is summarized 
in figure 14. Similarly we expect no parametric difference in the final equilibration times. 
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A Daughters neglecting re-splitting 

In most of section 4, when treating the split-off daughters of "hard" p ~ Q excitations, we 
made the simplifying approximation that q{9) is uniform in angle and nearly independent 
of an excitation's energy. We justified this approximation by arguing that we were really 
most interested in any new physics which comes to supplant g; nst as the dominant source 
of q. And in every case we found that the new q has at most order- 1 anisotropy. Therefore, 
at the moment when a new mechanism takes over controlling q, one is free to treat q as 
isotropic. But what happens before this moment? 

We will present at least a partial answer to this question, by looking at the evolution 
of particle occupancy and angular distribution, as a function of momentum p and time t 
for m < p < Q. We will treat tbroaden > t > 1/m, and we will also only try to understand 
p > k sp \it(t). We will only consider the case of an oblate momentum distribution. 

Further, we will neglect elastic scattering. We saw previously that this is a good ap- 
proximation for narrow angles and large momenta provided that g; ns t ~ 5~ 2 m 3 3> ^elastic ~ 
am 2 Q, which requires l + c + 3d > 0. For generic angles it requires qi QS t(0 ~ 1) ~ 5~ l m 3 2> 
am 2 Q, which requires 1 + c + d > 0. For the lowest-momentum excitations, p ~ m, it 
requires m 3 3> am 2 Q, which is true if 1 + c — d > 0. These restrictions, and the restriction 
p > k sp \\ t (t), must be taken into account when applying the results of this appendix. 

Our results are presented in figure 15. It is most convenient to write p = m5~ a and 
t = 6 /m, that is, to use m as the characteristic scale and powers of 5~ l to distinguish how 
much larger than m, 1/m the momentum p and time t are. In these units, Q ~ 5~ m with 
K = (1 — c + d)/{2d). Lines of equal K are lines of constant slope starting at c = 1, d = 0; 
the line c = 1 — 3d corresponds to Q ~ 5~ 2 m and the d = axis is K = oo. 

First, we determine the range of angles which split-daughters will take, as a function 

of p and t. The key facts we need are the following. Radiated daughters are naturally 

born in the narrow angular range of the hard parent population, and they then undergo 

transverse momentum diffusion to gain mean-squared momentum Ap 2 ~ qt. The range of 

angles will be 

if less than 1 

(A.l) 




The value of q is dependent on both p and 9 (see eq. (4.6) and eq. (4.8)): 

(777 \ 

— , m 2 p\ . (A.2) 

First consider excitations with p < m/5; then we always have q ~ m 2 p. At early times 
the angular range is 6 2 ~ qt/p 2 ~ m 2 t/p. These particles become isotropic on a time scale 

t ~ — r ~ n ~ — Tn . (A.o) 

q pm z m 
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Figure 15. The momentum-time plane for radiated daughters. Regions separated by solid lines 
have daughters experience the indicated values of q and fill the indicated angular ranges 9; in region 
1, q ~ m 2 p and 9 ~ mp~t2 . Regions separated by solid and dotted lines are numbered to indicate 
how the particle occupancy varies, with f(p,t) listed in eq. (A. 6). The figure is only valid before 
further splittings occur, that is, above the line k sp nt(t) (not shown, dependent on the ratio S/a). 
All lines have slope 0, \, 1, or —1 except the dotted line between regions 4,6 which has slope |. 



Next consider m/5 2 > p > m/5. Initially q 
when m 2 p ~ m s /59, which is at t ~ 5~ 2 p~ ] 



m 2 p and again 9 ~ mt 2 p 2 ; but this changes 
After this, q ~ m 3 /59, and we find 



qt 



m t 



r ^i _2 1 

3 rap si 3 



(A.4) 



-1 2 — '\ 
p m . 



p 2 p 2 59 

until 9 ~ 1 at the time scale t ~ 5~ 

Finally, if p > 5~ 2 m, then initially 9^5 the angular range inherited from the radiating 
parent. Enough momentum broadening has accumulated to overcome this value when 
eq. (A.4) returns 9^5, which is at t ~ 5~ 4 p 2 m~ 3 . After this the value of 9 is given by 
eq. (A.4) until t ~ 5~ l p 2 m~' i , when it is again order-1. These time and momentum scales, 
q values and angular ranges are summarized in figure 15. 

Regarding the occupancy as a function of p and t, we already saw that the formation 



time for a radiated daughter is t 
is tform ~ 5p2r77~2 for p > 5~ 2 m. 



2 

form 



p/q, which is tf Q 



~ 777 



for p < S 2 m and 



The rate of particle emission by a hard parent is 



T ~ a[l+f(Q)]/t{ orm (dtdp/p) and the density of hard parents is nhard 
powers of a cancel, so the number of soft excitations n p = dn/dlnp is 



777 3 Qt[l + /(Q)] 



p <c 5 777 , 



S^p-zmzQtll+fiQ)} p > <5 



777 , 



a 1 m 2 Q. The 



(A.5) 



Here [l+f(Q)] is a stimulation factor which is important when the hard particles are at 
large occupancy. 
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m 

1+/(Q) 



To convert these into an occupancy at the scale p, we use n p ~ p 3 6f(j>). The factor of 9 
is because not all of the angular range is filled. When the resulting occupancy is larger than 
Q[i.+f(Q)]/p, then joining (re-absorption) will compete with splitting and the occupancy 
will saturate. This occurs for t > p 2 m~ 3 for p < 5~ 2 m and for t > Sp^m~^ for p > 5~ 2 m. 
Therefore, in terms of the nine regions marked in figure 15, the typical occupancy will be: 

in p zQt? region 1 , 

m 3 p~ 3 Qt region 2 , 

Q/p region 3 , 

Q/p region 4 , 

( m 3 p~ 3 Qt region 5 , (A. 6) 

5~ l mzp~2Qt region 6 , 

53m 2 p~3Qt3 region 7 , 

o 3m2p 6 Qta region 8 , 

_ 5~ 2 mzp~2Qt region 9 . 
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